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This is an ebulliometer—being assembled by the master glass- 
blower on his own lathe at our Central Research Laboratory. 

The ebulliometer uses two vapor lift pumps such as those found 
in electric coffee percolators. One circulates a polymer solution 
while the other circulates boiling solvent—each liquid stream is 
directed over a temperature detector so sensitive that a difference 
of one-millionth of 1°C. may be sensed. This device allows the 
measurement of high molecular weights related to the physical 
properties of polymers. 

Such C-I-L research contributes to higher standards of per- 
formance in the plastics industry . . . to new and better plastic 
products . . . and—through better technical service—to greater 
success in the application of these new materials to modern living. 
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Laval University. 
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A Radioactive Tracer Technique for 


Particle Velocity Measurement 


in Solids‘Gas Systems’ 


W. H. GAUVIN?, I. 


To record the velocity history of a particle in 
complex problems of solid-gas flow with greater 
accuracy than is possible with existing methods, a 
technique was developed which in essence consists of 
following the flight of the particle under test—which 
has been previously tagged with a suitable gamma- 
emitting nuclide—by means of Geiger probes placed 
at predetermined points along the outside surface 
of the duct or reactor and parallel to its axis. Sharp 
timing pulses are thus obtained on the screen of a 
cathode-ray oscilloscope, which permit the estimation 
of the time elapsed between probe stations to less 
than 1/1000 second. Complete description of the 
electronic aspects of the method is given, and details 
of the particle preparation are outlined. 

A brief description is given of the application of 
the method to the determination of the rate of mass 
and heat transfer from particles in accelerated 
motion in a turbulent co-current gas stream and to 
the study of the effects of the turbulent flow para- 
meters on the drag coefficients of particles in highly 
turbulent fields. 





UNDAMENTAL studies involving the transfer of heat, 
F mass or momentum from particles entrained in a 
moving fluid share a common difficulty in that they 
require an accurate knowledge of the particle velocity 
history. In addition to being capable of good reproduci- 
bility, a suitable technique must allow for measurements 
to be taken simultaneously at many points in a flow field, 
which could be over a foot in diameter, without the intro- 
duction of any disturbance. Gasterstadt “') was able to 
determine the velocity of metal balls moving in a two-inch 
diameter glass tube by recording the change i in inductance 
which their passage w ‘ould pr oduce ina coil in contact with 
the circumference of the tube. This method, however, is 
limited to systems involving relatively large metallic 
particles being conveyed in relativ ely small conduits. 


Manuscript received January 6, 1959. 

2Department of Chemical Engineering, McGill University, Pulp and Paper 
Rescarch Institute of Canada, Montreal, Que. 

8Department of Chemical Engineering, McGill University, Montreal, Que. 
‘Department of Chemistry, McGill University, Montreal, Que. 
Contribution from the Department of Chemical Engineering, McGill 
University, Montreal, Que. 
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PASTERNAK*, L. B. TOROBIN?® and L. YAFFE*+ 


The triggering of a photoelectric circuit by particles 
has been commonly employed to obtain the settling 
velocities of particles in a fluid, but the need for an 
accurately focussed beam of light limits the method to 
single- -patticle studies, performed in transparent tubes of 
small diameter. An unsuccessful attempt to employ a 
photoelectric system to obtain the drag coefficients of 
freely-falling shapes i in air has been reported by Barker ‘) 
Ballistics experiments usually involve the perforation by 
the particle of a fine screen, a practice which would 
introduce large error in velocity studies at low Mach 
numbers. 

A method frequently employed in fluidized transport 
and pneumatic conveying experiments involves the 
instantaneous isolation of a section of the conduit by the 
simultaneous closing of two slide valves. The solids con- 
tent of the section is then related to the mean particle 
velocity in the system. Mitlin ‘*) has shown that measure- 
ments of this type are capable of great error, unless 
tedious precautions are taken. It would not be applicable 
to the measurement of velocities in the feeding, accelera- 
tion or curved sections of such a system. 

High speed cinematography is becoming increasingly 
porn with the advent of improved equipment. The 
latter is, however, costly and the analysis of results for 
non-steady motion is difficult since an individual motion 
camera, firing perhaps 100 ft. of film, is required for each 
experimental point. Khudiakow “) has shown that the 
method may not be used for single particles less than $-in. 
in diameter and at the same time in tubes larger than 2-in. 
The system has the further weakness of being limited to 
transparent conduits. A more extensive survey of particle 
measurement methods is given in a monogr aph on solids- 
gas flow “). 

The authors are at present engaged in a study of 
particulate behavior in turbulent fields. Two problems 
in particular require accurate knowledge of the particle 
velocity: the first is concerned with the determination of 
the rate of heat and mass transfer to individual shapes 
moving in accelerated motion through a turbulent co- 
current air stream, while the second involves a study of 
the effect of the turbulent flow parameters on the drag 
coefficients of particles in highly turbulent fields, in both 
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steady and accelerated motions. The methods of particle 
velocity measurement previously mentioned were all 
found to be unsuitable, and an accurate technique giving a 
relatively complete record of an individual particle history 
was therefore developed. 


Experimental method 


In essence, the method consists of following the flight 
of the particle under test—which has been previously 
tagged with a suitable gamma- -emitting nuclide—by means 
of Geiger probes placed at predetermined points along 
the outside surface of the duct and parallel to its axis. As 
the particle passes each probe it fires the Geiger circuit, 
producing a sharp timing pulse on the screen of a cathode- 
ray oscilloscope, to the vertical poles of which the probes 
are connected in parallel. The distances between pulses, 
suitably photographed, are referred to a signal produced 
simultaneously from a pulse generator, and permits the 
estimation of the time elapsed between probe stations to 
the nearest 1/1000 second. The results obtained are then 
plotted on a time-versus-distance graph or analyzed by 
numerical differentiation methods. 


Probe circuit 


In applying this method, a series of Vic-Tic end 
window Geiger probes, Model 631-56, equipped with 


IB85 counter tubes, were placed alongside the test duct. 
These were connected in parallel to the vertical poles of 
a Hewlett-Packard Model 130A cathode ray oscilloscope. 
The latter has a sensitivity of from 1 my./cm. to 50 v./em., 
and can be operated to sweep times as high as 15 sec./cem. 

The pulses from the Geiger circuit, together with the 
timing signals produced by a Wollensak model WF-311 
crystal- -calibrated pulse generator, were both recorded on 
Adox KB-14 35-mm. film as they appeared simultaneously 
on the screen of the nee A ty pical oscilloscope 
pattern is given in Figure The fine structure of the 
pattern, i.e. the fine pe: ak of e: . pulse, as well as the faint 
downward timing pulses, may not be easily discerned on 
the photograph as given, nevertheless, the film is easily 
read with the aid of a projecting microscope supplied 
with a calibrated stage. 


The complete equipment set-up required for the 
velocity measurements is shown diagrammatically — in 
Figure 2 

An Atomic Energy of Canada Limited Geiger and 
scintillation probe monitor, Model AEP 1903-S w hich had 
been modified by the authors, supplied the probes— con- 
nected in parallel with their required activation voltage. 
It is equipped w ith a counting rate meter which operates 
from an integrating circuit, indicating the amount of 





Figure 1—Typical oscillocope pattern. 
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Figure 2—Diagram of velocity measurement equipment, 


radiation intercepted by the probes per unit time. This 
rate increases as the particle approaches a probe, and 
becomes a maximum when the particle is directly in front 
of it. The characteristic pulse obtained for a given geom- 
etry and particle activation level will vary with the time 
constant of the integration circuit, and maximum of 
pulse definition is obtained by varying its capacitance 
value. Appropriate capacitors were inserted in the meter 
circuit to have time constants of 0.8 and 0.08 seconds, the 
latter being best suited for the particle activation levels 
employed. 

Since the circuit will also detect nuclear events 
resulting from cosmic radiation and ambient radioactive 
nuclides, it is necessary to remove this contribution by 
proper attenuation of ‘the pulse signal. This is accomp- 
lished by the inclusion of a variable capacitor into the 
range selection circuit of the monitor, and also by the 
a selection of the oscilloscope amplification level. 
Capacitors were added to the range circuit so that the 
instrument could register up to 100,000 counts per minute, 

The integrated signal, which had a full scale valve of 
40 mv., was taken from a point in the circuit immediately 
before the meter terminals which was 105 v. above ground 
level. It was passed through a variable voltage battery 
pack set to oppose the signal with 105 v., so that only the 
40 my. component carried through to the oscilloscope. 


The pulse generator provides a balanced output of 
approximately one v., and it is reduced to a level of the 
same order as the integration circuit pulses by means of 
the voltage divider illustrated in Figure 3. It will allow 
a calibration pulse height variation of from 0 to 10 mv. 
A 0.01 micro-farad capacitor has been placed in series 
with the oscilloscope input to eliminate parasitic back 
currents. 

The relatively inexpensive IB85 Geiger tubes were 
found to be quite adequate for the accuracy required it 
the authors’ investigational program. These have, aul 
ever, a very low efficiency for gamma ray counting 

(approximately 1%) so that the accuracy of the detection 
method can be further increased by the substitution of 
scintillation probes which have efficiencies of up to 100% 
These can be fitted with narrow slits and fed into a pulse 
height analyzer which would reject all pulses, except 
those which had passed through the slit. 

The timing technique described could be replaced by 
a series of decade counters. Each decade counter r. adjusted 
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Figure 3—Voltage divider. 


to start counting at the peak of a pulse, would operate 
between two probes, the first pulse opening the timing 
gate, and the second one closing it. The time interval 
would be indicated directly on ‘the counter panel, and 
the time intervals could be given to six figures if desired, 
provided that the particles were activated to a sufficiently 
high level. 

(b) Particle preparation 

The minimum particle surface activation level required 
by the equipment employed and experimental conditions 
was found to be approximately 20 mr./hr. In the case of 
celite particles down to |-in. in diameter, it was found 
that this material absorbed sufficient neutrons in a 72-hr. 
irradiation in the Chalk River N.R.X. reactor in a position 
having a. flux density of 8.6 x 10!2 neutrons/sec. x cm.” 
to bring it to this level. 

The natural neutron-absorbing ability of spheres hav- 
ing a wide range of densities and sizes, down to 500 
microns in diameter, cannot be relied upon to achieve the 
required activation level. These spheres were, however, 
readily prepared i in the following manner: Small irregular 
particles of the desired material ‘(celite, cork or pith have 
so far been used) were placed in a ball-making apparatus 
consisting simply of a circular race, two inches in diameter 
and a half an inch in depth and lined with abrasive strips 
of varying degrees of roughness. A tangential air jet 
issuing from a hypodermic needle w hirled the particles 
along the surface of the race, the attrition process 
shaping them into very accurate spheres in a period gener- 
ally not exceeding three minutes. The spheres were then 
soaked in a solution of iridium trichloride and dried 
under vacuum so that a fraction of a milligram of the 
salt remained in the interstices of the porous material. 
The particles were then sprayed with a concentrated 
solution of DuPont Zytel 61 Nylon resin, and dried at 
140 F. for three quarters of an hour. The particles were 
then restored to accurate sphericity by being re-run in 
the ball apparatus, operated without an abrasive strip. In 
this manner the radioactive tracer, when activated, is 


vd 


Figure 5—Velocity measuring equipment. 
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Figure 4—Probes in position along duct surface. 
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(c) Applications 

The accuracy and reproducibility of the method were 
intensively tested in connection with the two problems 
which have already been mentioned. The detailed findings 
will be discussed in subsequent articles. Only brief des- 
criptions of the experimental procedures will be presented 
at this time to illustrate typical applications of the velocity 
measurement method. 


AIR Vo= 56-0 ft/sec. 
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Figure 6 — Time 
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Figure 7—Initial behavior of a sphere falling in a stagnant 


fluid. 


To determine the rate of evaporation in accelerated 
motion from solid particles wetted with acetone, radio- 
active celite pe were dropped vertically from the top 
of a 14 in. pyrex gl ass line, 16-ft. in height, in the presence 
of a downward co-current turbulent flow of hot. air. 
Seven probes were installed against the outer surface of 
the insulation covering the glass line, two of which can 
be seen in position in Figure 4 , Which shows the lower 
portion of the mass transfer column, at the left. Not 
clearly seen is the receptacle containing water, at the 
bottom of the column, in which the falling sphere is 
caught, and the unevaporated acetone measured. Figure 
5 shows the electronic equipment with the oscilloscope in 
the back, its screen covered by black bellows to which 
the camera can be affixed. Shown in front, from left to 
right, are the pulse generator, the triggering circuit, the 
Geiger and scintillation probe monitor, containing the am- 
plification and integration circuit, and the v -ariable voltage 
battery pack. The ‘small black box on top of the oscillo- 
scope houses the voltage divider. Not shown is the trigger- 
ing device (located on top of the glass column) w hich 
starts the sweep of the oscilloscope simultaneously with the 
release of the particle. Typical of the time-versus-distance 
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data obtained is Figure 6, which describes the particle tra- 
jectories for a family of spheres, ranging in size from 0.15- 
in to 0.39 in., conveyed by air at 410°F. flowi ing at 56 
ft./sec . Interpretation of the complex mass transfer phe- 
nomena occurring in such a system could not be attempted 
if the contact times and continuous velocity history had 
not been determined. 


The second application which is currently being 
employ ed is concerned with the effect of high relative 
intensities of turbulence on the drag coefficients of spheres 
and shapes moving in a turbulent air stream. The spheres 
are fired singly in a manner which does not impart any 
rotation to them, through the grid inlet system and along 
the axis of an 8-in. i.d., 24-ft. high vertical turbulent wind 
tunnel which has been calibrated for the intensity, scale 
and spectrum values of the turbulent flow. In this manner, 
the effect of the turbulence parameters relative to the 
boundary layer and wake flow of the particle, as well as 
the effects of various rates of deceleration and acceleration 
are being studied. 

A particle dropping device has also been developed 
with which the drag on the sphere in the first few 
moments of flight has been studied, a typical time-versus- 
distance result of which is given in Figure ip 


All the results so far obtained follow smooth curves 
without any measurable deviations. No attenuation in 
sharpness was found when the system was made to 
correspond to a 12-in. diameter duct ‘surrounded by a 1-in. 
layer of thermal insulation. The probes can also be placed 
along any portion of a pipe bend or feeding section so 
that the average velocity of a multi-particle system can 
be obtained by the inclusion of a single particle in which 
a minute radioactive source has been implanted. 


It is hoped that the two typical illustrations which 
have just been discussed will suggest many more such 
applications for this simple method and thus help to 
advance the knowledge of the fluid mechanics of solids- 
gas systems. 


References 


Gasterstidt, J., Diss. Tech. Hochsch. Dresden (1922). 

Barker, D. H., Ph.D. Thesis, University of Utah (1951). 

Mitlin, A., Ph.D. Thesis, University of London (1954). 
Khudiakow, G. N., Invest. Akad. Navk S.S.S.R., Otdel, Tekh. Navk, 
7, 1022-1034 (1953). 

Torobin, L. B., and Gauvin, W. H., ‘‘Fundamental Properties of 
Solids-Gas Flow”, Pulp and Paper Research Institute of Canada 
Monograph (1956). 


( 
(< 
(+ 
( 


moe 


(5 


x * * 


The Canadian Journal of Chemical Engineering, June, 1959 





vel 
pu 
str 


pre 
Me 
rec 
the 


ty] 
fu 
ph 
iu! 
in 


ho 
SC] 


pr 
an 
pr 
be 
or 


te 





tra- 
15- 
56 
he- 
ted 


ad 


ing 
ive 
res 
res 
ny 
Ng 
nd 
ale 
er, 
he 

as 
on 


ed 
WwW 
Is- 


In 
to 


ed 
sO 
in 
th 





Principles and Operation 
of an Air Operated Mixer-Settler’ 


W. G. MATHERS? and E. E. WINTER? 


A mixer-settler for continuous countercurrent sol- 
vent extraction is described in which the mixing and 
pumping of the liquid phases are achieved by air 
streams. 

The operating principles of this mixer-settler are 
presented in theory and confirmed by experiment. 
Measurements include interface heights, capacity, air 
requirements and entrainment of aqueous phase in 
the solvent under various operating conditions. 

The main advantages claimed for this unit are. 

1) Mechanical simplicity and consequent free- 

dom from maintenance. 

2) Pumps and interface controllers are not re- 

quired between individual stages. 

3) Compactness for a given separation. 


1) Low energy requirements. 


HE nuclear power reactors to be built in Canada in 
Tike near future will be the thermal, heterogeneous 
type wtih a heavy water moderator and natural uranium 
fuel. The spent fuel from these reactors will contain 
plutonium, highly radioactive fission products and uran- 
ium. The plutonium is a valuable fuel and will be used 
in future power reactors, but it will not be used in the 
first Canadian reactors. Considerable work has been done, 
however, on processes and equipment to carry out the 
separations. 

At the present time the most economical and versatile 
process for the recovery and purification of the uranium 
and plutonium from spent fuel is a solvent extraction 
process. The highly radioactive nature of the solutions 
being handled in this process impose certain conditions 
on the extraction equipment to be used. 

(1) It should be essentially maintenance free. Main- 
tenance in a radio chemical plant is difficult and expensive. 

(2) It should occupy a relatively small volume for a 
given throughput and separation duty so as to reduce the 
amount of shielding. 


1Manuscript received October 30, 1958. 
2Atomic Energy of Canada Limited, Chalk River, Ont. 
Contribution from Atomic Energy of Canada Limited, Chalk River, Ont. 
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(3) A horizontal development is preferable to a 
vertical one because of the high cost of construction of 
shielded towers. 

At Chalk River emphasis has been placed on an air 
operated mixer settler known as the airlift contactor. 
This unit meets the three conditions given above. It is 
essentially a multi-stage, continuous, counter current, 
mixer settler in which the mixing and pumping of both 
phases from stage to stage are carried out by air streams 
introduced at the bottom of each mixer. Russell and 
Winter describe the development of an airlift contactor 
with a nominal uranium processing capacity of 8.5 tons 
per year “), The present work describes the design and 
operation of a three stage air-lift contactor with a nominal 
processing capacity of 100 tons of uranium per year. In 
an actual plant each contactor would have to have more 
than three stages but the data obtained from the three 
stage unit can easily be extrapolated to larger units. 

Although the air lift contactor was developed specific- 
ally for reactor fuel processing certain of its features 
might make it of more general interest. Besides those 
mentioned abov e, the power requirements, for the same 
conditions, are considerably less than those of a mechanic- 
ally stirred mixer-settler. This will be demonstrated 
further on in the paper. 


Contactor description 

\ drawing of the three stage contactor is presented in 
Figure 1. Referring to Figure 1 those chambers containing 
lift tubes are mixers and the others settlers. In operation, 
air streams are introduced at the bottom of each mixer 
and rise through the totally submerged lift tube. A circu- 
lation of liquid is thus set up w ithin the mixer and leads 
to a thorough mixing of the two phases. A schematic 
diagram of the flow through the contactor appears in 
Figure 2. In common with all mixer settlers the flow of 
the two phases is countercurrent between the stages and 
concurrent in any given stage. The two phases enter a 
given mixer from opposite directions, are contacted, and 
then flow into the settler through the mixed phase port. 
The phases are separated by grav ity in the settler and the 
light phase proceeds to the next mixer through a baffled 
rectangular weir. The function of the baffle is to minimize 
a backward flow of mixed phases into the settler by 
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Figure 1—Three stage air lift contactor isometric view. 


splashing. The heavy phase proceeds in the opposite direc- 
tion through an underflow line into the next mixer. 
Referring to Figure 1 it can be seen that this line rises into 
the mixer. It will be demonstrated that the height of 
this riser serves to determine the interface position for 
a given set of operating conditions. 

The experimental contactor was constructed from 
stainless steel with plates of the transparent plastic Hom- 
alite on the front, top and back. Thus visual observation 
of the contactor operation was possible. 


Principles of operation 


In a mixer settler with inter stage connections as 
shown in Figure 2, all liquid levels and interface heights 
can be determined in terms of known quantities provided 
the an requirements are met: 

(1) The mixing is “perfect” , Le., the ratio of the two 

ities in a mixer is everywhere the same. 

(2) Pressure drops due to flow throughout the appa- 

ratus are negligibly small in comparison with 
various static pressures. 


LIGHT PHASE 
OuT 








~ 





SETTLER 
' 


HEAVY PHASE 
IN 


SETTLER 
2 


(3) The dynamic pressure w hich is due to the mixin 
circulation at the point where the heavy phase 
enters the mixer and at the mixed phase port must 
not differ appreciably from the static pressure at 
those points. 


Criterion (2) above is a question of proper design, 
Whether or not criteria (1) and (3) are satisfied for a 
particular mixer settler can only be settled by test. 

Referring to Figure 2 the total depth of liquid in 
settler 1 will be uniquely determined by weir geometry 
and light phase nei rate. The difference in lev a between 
settler 1 and mixer Ah,, can then be determined by an 
hydraulic balance aan the mixed phase port. If the 
mixing is “perfect” the ‘fluid in the mixer will have an 
average density given by 


Pav, 1 = fu,1 puoi + (1 — fe i)pr, 1 


where f,; is the volume fraction of heavy phase in the 
mixer. 

p is the fluid density. 
Using this expression the balance through the port gives 


Ah; = Ee acee hy = JN yen nes (1) 


1—f 
1+ omen i 


fi, 1 


where Ah is the difference in the liquid height between 
the settler and mixer h,,h, refer to heights as shown in 
Figure 2. 

Examination of Equation (1) shows that Ah will 
always be positive, i.e., in a given stage the total liquid 
depth in the settler will always be greater than that in 
the mixer. If Ah, is large enough such that the liquid 
level in mixer 1 lies below the crest of the weir connecting 
settler 2 and mixer | then the flow through this weir will 
be unimpeded and the total liquid depth in settler 2 will 
be the same as that in settler 1 since the weir geometries 
and light phase flow rates are the same for both stages. 
This argument can be extended to show that, in a multi- 
stage contactor provided individual Ah’s are large enough, 
the total liquid depth in all settlers will be the same. 
Since there is no hydraulic gradient in the settlers it can 
be seen that the mixing, in effect, pumps the light phase 
through the contactor. Requiring the liquid depth in a 
given mixer to be below the crest of the weir emptying 
into it is equivalent to saying that Ah must be greater 
than the height of liquid over the weir crest in the 
settler. Ideally then one could use Equation (1) and the 
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Figure 2— Schematic flow diagram for the three stage 
mixer settler. 
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app opriate weir formula to obtain a crest width which 
would prevent the appearance of an hydraulic gradient 
for any given set of oper ating conditions. In practice, 
how ever, the term f,,; in Equation (1) may not be known 
(as will be shown later) and weir formulae breakdown 
for the small depths considered here so that a considerable 
safecy factor should be incorporated into any such pro- 
cedure. There is nothing wrong, for instance, with having 
weirs which are wider than they need be. 


it should be mentioned here that Equation (1) will 
only apply when the interface in each settler lies oe 
the mixed phase port between the mixer and settler. | 
it should be above this port the liquid depth in the a 
can become greater than that in the settler thus impeding 
the light phase flow into the mixer. This condition can 
set up an hydraulic gradient in the contactor. It will be 
shown below under what condition this can take place 
in the model under consideration. 

To determine the interface height in settler 1 and hy- 
draulic balance is made from settler 1 through the heavy 
phase v underflow line to mixer 2 and through the mixed 
phase port to settler 2. This gives 


hi, 1 pu 1+ (hy — hi 1) pr 1 = hp Pu, 1 + 
(h, — hr) (fu, 2 pu, 2 + (1 — fi.2) pe, 2) + (he — he) pr, 2 (2) 
where 


h’s refer to heights as shown in Figure 2. 
Equation (2) can be rearranged to solve for the inter- 
face height and gives: 


Pu, 2 — PL, 2 


hi,1 = hofu, 2 + 
me 2. — Be. 1 
" pu, 1 — (fu, 2 pu, 2 + (1 — fu, 2) pr, 2) 4 
Bn; 2 By. 
h, rTM mee Le Wie ete (3) 
ni. By. 


Equation (3) is a general expression for interface 
height for the case where the density of both phases 
is varying from stage to stage. Some special cases can 
now be considered. 

When there is no density change between stages in 
either phase Equation (3) becomes 


hi. 1 = fu, oh, a (1 — fu » hr nt 3 (4) 


and, in particular, if the stages are connected such that 
h, = O then 


hi, ge Bales, oe ‘ : (5) 


fu, 2 Will always be between O and 1 so that, in this case, 
h, > h, , > O,i.e., the interface lies between the required 
limits. 

In the experimental work that follows the validity : 
the assumptions leading to Equation (3) were tested i 
a system where no density changes took place. aunt 
(+) should then apply. 

On examination of Equation (3) it can be seen that 
the first term on the right hand side will always be 
positive and the second term will also be positive except 
for the unusual case where the heavy phase in stage | 
is less dense than the mixed phases in mixer 2. It is very 
unlikely that this situation would ever be met in practice 
so the second term on the right hand side of E quation (3) 
can also be assumed to be alw ays positive. The sign of 
the third term can be either positive or negative depending 
on whether the light phase density increases or decreases 
in the direction of its flow. If the light phase density 
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increases in the direction of light phase flow the third 
term becomes negative and the interface is lower than 
it would be in the absence of density change. The inter- 
face can be maintained above the heavy phase underflow, 
however, by a suitable choice of hx. 

A more serious condition can exist if the light phase 
density decreases in the direction of its flow. Here the 
interfaces will be high and in some systems could rise 
above the mixed phase port and cause an hydraulic 
gradient to appear. When this situation exists h, would 
be made equal to zero in order to keep the interfaces as 
low as possible and by requiring that h, , < h, the follow- 
ing inequality is obtained from E quation (3). 


h, J 
+ — G1 < Pu, 1 — Pr, 1 — fu, 2 (Pu, 2 — Pr, 2) 


h, 

This inequality provides a criterion which must be 
met if the stages in a multi stage contactor are to be 
connected as in Figure 2. 

In the above material the treatment was carried out 
for the first and second stages of the contactor but it will 
be realized that it is just as valid for the nt’ and n + 1th 
stages. 

The interface height in the settler where the heavy 
phase leaves the contactor must be controlled externally. 
In this work an adjustable height weir was used. 


Interface level 


The validity of the assumptions leading to Equation 
(3) in which interface height within the contactor is 
predicted was subjected to test. The heavy phase fed to 
the contactor was 6% NH,NO, and the light phase was 
a 20% V/V solution of tributyl phosphate in Soltrol* 
(20% TBP-—Soltrol). There is no density change in this 
system so that Equation (4), should apply. The interface 
h; was measured in settler 2 and the phase ratio in mixer 
3 was determined by drawing samples of the mixed phases 
through sample taps. The results are summarized in Table 
h. ‘ 

TABLE 1 
COMPARISON BETWEEN OBSERVED AND CALCULATED 
INTERFACE HEIGHTS 
he = 4.75 ins., h, = 11-5 ins. 


Heavy Phase Interface Height 


Heavy Phase Light Phase Volume Fraction 


Flow Rate | Flow Rate in Mixer ths tte) 
(1./min.) (1./min.) (Measured ) 
fu Obs. Cale. 

0 1 | 1.0 0.100 ao 5 43 ; 
0.5 1.0 0.265 6.5 6.53 
1.0 1.0 0.365 eae esa 
2.6 1.0 0.496 8.25 8.11 
3.0 1.0 0.590 9.44 8.75 
a5 1.0 0.740 10.62 9.74 


This data is in good agreement with the theory except 
for the two final runs where, at an aqueous flow rate of 
3 liters/min. the pressure drop through the underflow 
line has started to become appreciable. 

An important result from the data of Table 1 is the 
fact that the measured fraction of heavy phase in the 
mixer departs considerably from the heavy phase fraction 
of the total flow fed to the mixer. Thus in the second last 


®Soltrol is the trade name of kerosene fraction. 
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Figure 3—Schematie diagram of the flow pattern between 
the mixer and settler. 


run although 75° of the total feed to the mixer was 
heavy phase only 59° of the mixed phases in the mixer 
was heavy phase. If one assumes “perfect” mixing, i.e., 
no point to point changes in phase ratio within the mixer 
and a positive flow of the mixed phases through the port 
into the settler then simple material balance considerations 
show that at the steady state the volume fraction of heavy 
phase in the mixer should equal volume fraction of heavy 
phase fed to the mixer. 


The assumption of “perfect” mixing was checked by 
drawing samples of mixed phases from five different 
points spaced at intervals between the top and the bottom 
of the mixer. All five samples showed identical phase 
ratios so that the assumption of “perfect” mixing is sound. 
One is therefore forced to the conclusion that there 
existed a reverse flow of light phase through the upper 
part of the port connecting the mixer and the settler. 
This flow pattern is illustrated schematically in Figure 3. 
If a reverse flow of light phase existed as shown it should 
be possible to eliminate it almost completely by position- 
ing a baffle where the dashed line appears in Figure 3. 
This was done while the contactor was in operation and it 
was observed that the fraction of heavy phase i in the mixer 
increased until it became the same as the fraction fed thus 
proving the existence of a reverse flow through the un- 
baffled port. From the design standpoint, then, the term 
fy, in Equations (1) and (3) will only be known as an 
upper limit unless the mixed cam port is baffled. 


This recirculation of light phase will have the effect 
of making the light phase retention time in the mixer 
somewhat greater than that of the heavy phase. 


Contactor effectiveness 


It has been mentioned above that a solvent extractor 
for a radiochemical process should occupy a relatively 
small volume for a given extraction duty and throughput. 
A relatively low hold up of radioactive solution will 
serve to reduce both the extent and the thickness of the 
shielding necessary. 


Coplan, Davidson and Zebroski‘*) suggest that a 
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measure of extractor compactness which they call con- 
tactor effectiveness is given by 


r NF 
Ey = — 
\ 
where E, is the contactor effectiveness 
N is the number of theoretical stages for a given 
system 
Fis the total flow rate 
Vis the contactor volume 
E, has the units of reciprocal time and is a direct 
measure of the relative compactness of a given extractor. 
In order to assess the effectiveness of the airlift con- 
tactor a run was made in which acetic acid was extracted 
into hexone (methy] isobutyl ketone). This system was 
used by Coplan et al to test the KAPL “pump- mix” mixer 
settler (2) and by Chantry and Berg in a packed column 
with and without pulsing (), Thus a comparison was 
possible with these three types of extractor. 


The conditions of the run are presented in Table 
In this run the mixer volumes were reduced from 10 
liters to 5 liters, using inserts, in order to obtain the same 
mixer retention time as was used by Coplan et al (2), 


TABLE 2 

EXTRACTION OF ACETIC AcID INTO HE XONE 
Aqueous feed 
Solvent feed 
Aqueous feed rate 
Solvent feed rate 
Air rate/stage 
Mixer volume 
Mixer retention time 
Settler volume 


- 0.969 acetic ac “id cole with hese »*xone 
-hexone saturated with water 
- 3.48 liters/min. 
6.90 liters/min. 
- 0.30 cfm 
- 5.0 liters 
- 0.48 mins. 
- 10.0 liters 


The results of the run are presented as a McCabe- 
Thiele diagram in Figure 4. The equilibrium curve was 
obtained by equilibrating s samples of the feeds in various 
volume ratios. The operating line was drawn through the 


lowest point with a slope “equal to aqueous to solvent 
volume feed ratio. 


The data yields an average Murphree stage efficiency 
(based on aqueous phase analy ses) of 93% . Coplan et al 
report an efficiency of 95 + 5 for te mechanically stirred 
“pump-mix” settler with the same mixer retention time 
and phase ratio. Thus the degree of mixing was com- 
parable in both units. 

In Table 3, the effectiveness of the airlift contactor in 
this run is compared with three other types of extractor. 


TABLE 3 


COMPARATIVE EFFECTIVENESS OF VARIOUS EXTRACTORS 


Contactor 
Effectiveness, 


2+ 

: Benterser . ee (hours!) 
PAGO AGA ES) so o.5 bose Rise NW wears i 
Pulsed Packed Column (3) ........... 33 
Pump-Mix Mixer Settler (2)............ 27 
Fat (Ait Content oii nike sw cuneate. 39 


°A higher value means a smaller volume is required for a given duty. 


Too much significance should not be attached to the 
numerical values in Table 3 since they were arrived 
for particular runs and do not necessarily represent maxi- 
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Figure 4—Extraction of acetic acid from water into hexone. 


mum values for the various units. One can say, however, 
that pulsing a packed column will reduce the volume 
necessary for a given duty by a factor of about three 
and that the pulsed packed column, the “pump- mix” 
mixer settler and the airlift contactor are close to equiv- 
alent in their volume requirements. 


Air requirements 


An objection to the use of the airlift contactor for the 
processing of radioactive material is that the air stream 
leaving the contactors will necessarily contain small 
amounts of entrained radioactive liquid. If this air is to 
be discharged to the atmosphere it must first be decon- 
taminated to an acceptable level. It is obvious that this 
decontamination problem will becomes less severe at 
lower operating air rates. Consequently it is of interest 
to know, for a given set of conditions, the minimum air 
rate which will give satisfactory physical operation of 
the contactor. The contactor is considered to be operating 
satisfactorily when the mixers are filled with mixed phases 
and there exist no point to point changes in phase ratio. 
This condition was quite easily identified by visual obser- 
vation. 


Some measurements of minimum air flow rates per 
stage are presented in Figure 5. Water and 20% TBP- 
Soltrol were fed at a ratio of unity and at varying rates, 
as shown. It is a peculiarity of the airlift contactor that 
the minimum air rate required to establish satisfactory 
operation after an interruption in the air supply has 
allowed the phases in the mixer to settle (referred to in 


Figure 5 as a minimum start up air) is less than the air, 


required to maintain proper operation (minimum oper- 
ating air). The reason for this behaviour is that, on start 
up, ‘heavy y phase only is being pumped through the lift 
tube. Initially, then the pumping t takes place against a 
certain head. Once proper operation has been established 
the liquid in the mixer becomes essentially homogeneous 
as regards density and the pumping head disappears. As 
a consequence the pumping rate through the lift tube 
will increase to a value greater than that necessary to 
maintain good mixing and the air rate can then be reduced 
without affecting the mixer operation. 


Two sets of curves are presented in Figure 5 showing 
minimum start up of and operating air rates per stage 
both with the lift tubes in place and with them removed. 
It can be seen that, while the air required to maintain 
good mixing is nearly the same in both cases, less air is 
required to start up when the lift tube is removed. Under 
these conditions, then, there is some advantage in not 
using a lift tube. However, when a lift tube is used, a 
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Figure 5—Air flow requirements per stage for various 
conditions. 


more ordered circulation of liquid within the mixer 
results. In certain cases this is essential for proper mixer 
operation. For instance, it has been observed here that, 
although the 6” x 4” x 26” mixers of the airlift contactor 
could be operated satisfactorily over a wide range of 
conditions without lift tubes, when the dimensions were 
reduced to 6” x 2” x 26”, using inserts, the contactor 
could only be operated with the lift tubes in place. 


From the measurements of Figure 5 the magnitude of 
the air decontamination problem to be expected in a 
radiochemical process can be evaluated. For instance it 
can be estimated that, in a plant handling 100 tons of 
irradiated uranium per year by the TBP process, approxi- 
mately 2 cubic feet of air per minute will be contacted 
with highly radio-active liquid. Pilot plant experience in 
this laboratory has shown that there is no difficulty in 
decontaminating an air flow of this order to a level well 
below the tolerable limit ‘*. 


Energy requirements 


If the inlet air pressure is known for the data of Figure 
5 the energy delivered to the liquid in the mixer can be 
obtained. This was measured and found to be essentially 
constant at 1.6 in. Hg. gauge over the range of the data. 
The energy delivered by this pressure head (the velocity 
head is negligibly small) in horsepow er can then be 
obtained by multiplying the air rate in cfm by 3.4 x 10°. 


Referring to Table 2 it can be seen that, in this run, 
the energy input was 10-* horsepower per stage or 9.6 x 
10° horsepower per stage per liter/min. of total flow. 
Overcashier, Kingsley and Olney give a value of 1.5 x 
10-8 horsepow er per liter/min. of total flow for the same 
mixer retention time and efficiency in a mechanically 
stirred mixer (5), Thus the air lift contactor requires less 
energy by a factor of 15 for this case. 


Aqueous entrainment in the solvent phase 

A rough measure of the effect of air flow rate on the 
interfacial area per unit volume can be obtained by 
measuring the aqueous entrainment in the solvent stream 
at varying air rates. This experiment was carried out and 
the results are presented in Figure 6. The lowest air rate 
used corresponded to the minimum operating air rate for 
those conditions. It can be seen that, over the range of 
the measurements, the amount of entrainment does not 
depend very greatly on air rate. A possible interpretation 
of this data is that although an increasing air rate increases 
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Figure 6—Effect of air rate on entrainment. 


the rate of circulation of liquid within the mixer it does 
not affect the degree of dispersion of the two phases to 
any great extent. 

Some measurements of aqueous entrainment in the 
solvent phase as a function of throughput at constant air 
rate are presented in Figure In the lower curve two 
sections of York mesh packing were placed in each 
settler, one above the mixed phase port and one below. 
In the upper curve no packing was used. 

In steady operation a band of unsettled emulsion was 
observed at the interface. The sharp increase in entrain- 
ment at the right hand end of these curves occurred 
when the top of this band rose to the light phase weirs 
and overflowed. It can be seen that, with packing, the 
settler capacity before this point was reached was in- 
creased by a factor of three. Also, for a given throughput, 
the use of packing decreased the amount of entrainment. 
Conclusions 

The airlift contactor meets the requirements of a 
solvent extractor for radiochemical processing. Viz: 

(1) It is relatively compact and is developed horizon- 

tally rather than vertically, 


It has no moving parts and consequently requires 
little or no maintenance. 


(3) A multistage unit does not require pumps or 

interface controllers between the stages. 

In addition to the above features the power require- 
ments for a given duty are low in comparison to a 
mechanically stirred mixer settler. It is also a versatile 
unit in that it can be operated over a wide range of 
throughputs and phase ratios. 
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Figure 7—Effect of flow rate on entrainment. 


The chief objection to the use of the airlift contactor 
in a radiochemical process is the contamination of the 
air with entrained radioactive liquid. Experience has 
shown, however, that it is a relatively easy matter to 
decontaminate the air to a level well below the tolerable 
limit. 


Nomenclature 


f is the volume fraction of one of the phases in the mixer 

h’s are heights as given in Figure 2. 

Ah is the difference in liquid height between the settler and the 
mixer 

pis the fluid density 


Subscripts: 


1, 2, 3 refer to stages 1, 2, 3 respectively 
H,L_ refer to the heavy and light phases respectively 
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Mixing Efhciency Determinations 
for Continuous Flow Systems’ 


A, CHOLETTE? and LEONCE CLOUTIER? 


Models are defined for various mixing conditions, 
in continuous flow systems. Differential equations are 
derived which take into account an effective volume 
of mixing, possible short-circuiting, hold-up time of 
the system, partial displacement or piston flow. 


The values of the different factors contained in 
the integrated equations can be determined experi- 
mentally by the particular response of a given system 
to a sudden change in composition of the feed. 

A correlation of the effective volume of mixing 
and the agitator r.p.m. is presented. 


A design procedure is suggested for tank flow 
reactors when partial mixing occurs. 


IxING, which has been referred to as one of the most 
M intractable of all the unit operations of chemical 
engineering (1), has been the object of more and more 
fundamental work over the past years, partly perhaps on 
account of the recognized complexity of the subject (?). 
Many of the investigations in this field have dealt with 
batch systems and with the power requirements for 
mixing (3,4. 5,6), The ever increasing conversion to con- 
tinuous-flow processes in the industry has introduced 
various other aspects to the problem for others to = 
Among the numerous avenues of approach followed, 
few can be found dealing with the efficiency of hie 
In spite of all the work done in this particular field, much 
more remains to be accomplished as indicated in recent 
chemical engineering works ‘!:*,7), 

The present investigation was undertaken with the 
purpose of providing a theoretical basis for evaluating 
the performance of mixers. The method developed for 
measuring mixing efficiencies applies to continuous-flow 
systems but it is intended to extend it to batch systems as 
well. After presenting a quantitative treatment of various 
factors likely to be encountered in actual mixing, an 
experimental study relating to part of the theoretical 
investigation is exposed. 


Previous work 


MacMullin and Weber ‘8) were among the first to 
present an analytical study of the behavior of perfectly 
mixed systems, for continuous-flow mixing vessels in series. 
They, and others (% 1), used the term “short-circuit” or 
“mathematical by-passing” to denote conditions involving 
probability considerations, applicable to well agitated 
tanks in which the contents are substantially homo- 


1Manuscript received January 2, 1959. 

2Department of Chemical Engineering, Laval University, Quebec, Que. 
Contribution from the Department of Chemical Engineering, Laval Uni- 
versity, Quebec, Que. 

Based on a paper presented at the Joint A.I.Ch.E.-C.I.C. Chemical 
Engineering Conference, Montreal, Que., April 20-23, 1958. 
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Les auteurs définissent des modéles pour diffé- 
rentes conditions d’agitation dans les systémes a 
écoulement continu. Ils déduisent des équations diffé- 
rentielles qui tiennent compte d’un volume efficace 
d’agitation, de court-circuits possibles, du temps de 
rétention du systéme et de l’écoulement par déplace- 
ment ou “frontal”. 

Les valeurs des différents facteurs contenus dans 
les équations intégrées peuvent étre obtenues 4 partir 
de la réponse particuliére d’un systéme donné a un 
changement subit de la composition de l’alimentation. 


Les auteurs présentent une corrélation entre le 
volume efficace d’agitation et la vitesse de l’agitateur. 

Ils proposent une méthode de calcul pour les 
réacteurs en régime continu lorsque l’agitation est 
partielle. 


geneous, whereby some of the molecules in the feed 
stream appear almost immediately in the effluent as a 
result of good agitation. This may be confusing because, 
as pointed out by Colburn “), the term “short-circuit” 
should refer rather to the flow pattern in a tank where 
some of the feed is channeled to the outlet, reaching the 
latter without being mixed. It later occurred to Mac- 
Mullin “2) that short-cut or channelized flow between 
inlet and outlet would indeed alter conditions in a tank. 
He suggested, rather qualitatively, that it could perhaps 
be possible to evaluate the efficiency of a stirrer by 
studying the conversion rate of a reacting system whose 
kinetics are known. 

Statistical methods have been suggested to determine 
blending efficiencies ‘') for batches of varying concen- 
tration fed into a stirred tank. 


Danckwerts ‘!#) has introduced the concepts of ‘hold- 
back’ and ‘segregation’ in order to compare actual to ideal 
mixing conditions in continuous flow systems. The con- 
cept of hold-back relates to the residence time of various 
elements in the system; the author cautions: “the whole 
F-diagram (and possibly other information) is needed if 
the actual performance of the reactor is to be calculated”. 
The notion of real ‘short-circuit’ is mentioned but none 
of the F diagrams shown take it into account. As to the 
concept of ‘segregation’, also obtained from F-diagrams, 
it is introduced to give an indication of the efficiency 
of mixing in a system. While it points: to the degree of 
departure of a system from perfect mixing, it does not 
indicate fully how the system actually behaves. The H 
and S values thus defined involve rather tedious deter- 
minations besides being questionable as to their real 
advantage. 
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Figure 1—Perfect mixing: the feed is dispersed instan- 
taneously. The composition of the outlet stream is the 
same as that of the tank contents. 


Factors involved 

With two miscible liquids, in a batch process, a 
perfectly homogeneous mixture can be obtained when 
agitation, with the help of diffusion if necessary, is main- 
tained for a long enough period of time. In a continuous 
flow system, howev er, it may be that perfect mixing can 
never be attained. There can be stagnant regions in the 
reservoir, part of the feed may go directly to the outlet, 
piston or plug flow may take place. In order to study the 
behavior of such a system and the efficiency of mixing, 
one can consider the tank or reservoir to act as a damper 
to variations in feed composition. The mathematical 
treatment of this problem thus becomes similar to one 
performed on a capacity in a servo-mechanism circuit, 
where a transfer function exists, from which can be 
obtained transient or permanent responses. 

Three factors that may be met in actual mixing con- 
ditions are taken into account. 


(a) Effective volume 

For a given agitator, part only of the total tank 
contents may be well mixed the rest being taken up by 
stationary or stagnant zones. That part of the total 
volume where mixing is complete will be referred to 
as the ‘Effective Volume of Mixing’: it is the part which 
acts as a damper to any variation in feed composition 
and which is considered in determining the transfer 
function. A somewhat similar concept has been used by 
Germain“) for evaluating contact time in a tubular 
reactor: the author replaces the reactor volume, in which 
the — is variable and passes through a maximum 
of T,,, by a smaller ‘equivalent volume’ in which the 
temperature is uniform, at T,,. 

MacDonald and Piret °) have observed qualitatively 
that dead volume can retain solution for indeterminate 
periods of time within a vessel. Danckwerts ‘'4) also 
considered dead space in a mixing tank when introducing 
the concept of hold-back. 


(b) Short-circuit 

Depending on the relative position of the inlet with 
respect to the outlet, part of the feed may go directly 
to the latter without being mixed, causing channeling or 
a short-circuit. Only the remaining portion is thus avail- 
able for mixing. This situation is unlike that arising from 
the definitions of short-circuit introduced by some 
authors as mentioned already. 

MacDonald and Piret '*) report that they have ob- 
served channeling, in the course of experimental deter- 
minations, especially at low levels of agitation. 


(c) Piston flow 


In many systems, depending on conditions, part of the 
4 
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Figure 2—Partial mixing and short-circuit: a fraction of 
the total volume is perfectly mixed; part of the feed goes 
directly to the outlet. 


feed, or all of it, may give rise to piston or plug flow. 
Inasmuch as the liquids i in the stagnant and well agitated 
zones are not separated by a membrane, fluctuations or 
circulation may arise. Liquid in the stagnant zone may 
be replaced by some from the well agitated zone or from 
part of the feed. Such a circulation may go on steadily 
in the form of piston or plug flow. 

Nagata and others ‘!7 18.19) have studied mixing in a 
tubular reactor where piston flow exists, but in the 
absence of any short-circuit. 

Mixing problems met in actual practice may involve 
different combinations of these three factors: effective 
volume, short-circuit, piston flow. Some particular com- 
binations, more likely perhaps to be encountered than 
others, will now be considered. 


Theory 


As for a servomechanism element, the behavior of a 
mixing system can be studied by means of the answer to 
either a step, a linear or a harmonic demand. One only 
of the answers being necessary to determine a transfer 
function, the first or step demand has been chosen on 
technical grounds suggested by the experimental deter- 
minations to be made. For the different systems con- 
sidered, the following are the answers to a negative step 
demand. 


1. Perfect mixing 


The contents of a reservoir, Figure 1, occupy a volume 
of V cubic feet; the initial concentration is c, lbs./ft.*. 
q cu. ft./hr., of composition cy, are assumed to be fed 
continuously into the reservoir. The feed composition cy 
is changed suddenly at time t = O from c, to O, creating 
a step demand. At any given time t, mixing being perfect, 
the tank contents are uniform throughout and of com- 
position c, which is also that of the effluent. 


A material balance for a time element dt gives: 
qcrdt = qcedt + d(Vc) 


The feed composition, cy, being then zero, and the 
volume of the liquid in the reservoir being constant, 
Equation (1) can be written: 


de ‘ q O 
c= 

dt V 

Upon integration one obtains: 


q 
- t 
V 
/ 
c/Co = 
When plotted on semi-log paper, as in Figure 1, this 
equation gives a straight line. 
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Figure 3—Partial mixing and piston flow: a fraction of 
the total volume is perfectly agitated; piston flow takes 
place before or after mixing, as shown. 


2. Partial mixing and short-circuit 

As shown in Figure 2, only a fraction m of the total 
volume is assumed to be well agitated. A fraction n of the 
feed enters the zone of perfect mixing while the other frac- 
tion (1 —n) is assumed to be short-circuited directly to the 
outlet. The concentration, at time t, in the zone of perfect 
mixing is assumed to be c’. The material issuing from 
this zone is mixed with that portion of the feed which 
short-circuits the system, the mixture of the two giving 
a liquid of composition c at the outlet. One thus obtains: 


nqcrdt = nqe’dt + mVde’ 
or, 
de’ nq 
-c’ =O 
dt m\ 
nq 
mV 
were Me See cowneeawes (2) 
But, 
(1 — n) qcrdt + nqe’dt = qcedt 
and, 


cr = O.). c = ne’ 


Equation (2) then becomes: 


nq 
= t 


mV 
c/Co = ne 
Taking natural logarithms: 
nq 
Inc/cg = inn — ie 
m\ 


Thus, plotting c/c, vs qt/V on semi-log paper as in 
Figure 2, the value of n can be read directly on the con- 
centration ratio axis, at t = O. The straight line has a 
slope of — n/m, from which can be derived the value of 
m, once n has been determined. 


3. Partial mixing and piston flow 


With a fraction m of the total volume well agitated, 
the feed can be considered to give piston flow in two 
different ways, as shown in Figure 3: 

(a) The feed is introduced into the zone of perfect 
mixing, displacing an equal volume into the stagnant 
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zone: this liquid, in turn, displaces the other already 
present, by the mechanism of piston or plug flow. 

(b) The feed, introduced into the stagnant zone, dis- 
places an equal volume from the latter into the zone of 
perfect mixing which, in turn, supplies an equivalent 
volume of liquid to the outlet. 

In either case, the composition of the iiquid at the 
outlet would remain constant at c, for a period of time 

(1 — m) V/q, which represents the time required 
to displace all the material out of the stagnant zone at the 
prevalent rate of flow. The variation of the outlet com- 
position with time, in either case, is obtained in the 
following way: 


qcrdt = qe’dt + mVdc’ 


de’ c 
— . c’ =O 
dt m\ 
aet 
mV 
c’'/cp = e 
: . (1—m)V 
But, c = c’ at time: t — oe 
o q 
1 qt 7} 
= [3 al ah 
-). C/Op = € 
, 1 | qt E 
Inc/cg = —- — —~ — (1 — m) 
m]| \ 
Plotting this equation on semi-log paper gives a 
straight line of slope — 1/m, as shown in Figure 3. The 


value of m can also be obtained from the intercept at 
c/é, =}. 


4. Partial mixing with piston flow and short-circuit 

Among the several possible combinations, the follow- 
ing three have been considered. 

(A) In addition to having a fraction m of the total 
volume well mixed and a fraction n of the feed going 
to the zone of perfect mixing, a fraction p of the feed 
gets into the stagnant zone from where it displaces an 
equivalent volume by piston flow into the zone of perfect 
mixing, as shown in Figure 4. The outgoing liquid, of 
composition c at time t, is made up in part of liquid of 
composition c’ from the zone of perfect mixing, and of 
the fraction (1 — n — p) of the feed which short-circuits 
the tank. Two periods have to be considered, before and 
after time t = (1 — m) V/p q, during which all the 
liquid originally in the stagnant zone is displaced. 


a) O<t <(1 — m) V/pq 





Figure 4—Partial mixing with piston flow and_ short- 
circuit: part of the feed displaces liquid from the stagnant 
region into the zone of perfect mixing. 
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Figure 5—Partial mixing with piston flow and short-circuit: 
part of the feed displaces liquid from the stagnant region, 
around the zone of perfect mixing. 


During this period, material of composition cy = c, 
enters the zone of perfect mixing from the stagnant zone. 
The following relations are obtained: 

nqcrdt + pqe.dt = (p +n) qe’dt + mVdc’ 


de’ (p+n)q , pac 
1; : Cc = = 
dt m\ m\ 
— Ptnae 
v 
x p , m 
c’/c, = + Ke 
p+ nh 
Since 
‘ n 
c’=cqatt=O, Ke= 
pt+no 
_ (p+n) ate 
A p n mV 
c’/Co = + e 
pt+no pt+no 
c=(p+n)c’ 
(p+n)q 
mV 
.). C/Co = p + ne 


b) t > (1 — m) V/pq 


All the material originally in the stagnant zone has 
now left and liquid of composition cy = O enters the 
zone of perfect mixing from the stagnant one. 


‘ Ic’ 
ae ip +B) 4). 
+ - c’ =O 
dt m\ 
(p+n)q 
mV 
c’ = Kie 
(p+n)q 
- : t 
mV 
c/co = Kye 
‘ (1—m)V 
Now, at time: t = 
. Pq 
(p +n) (1 m) - Oa n) (1 — m) 
mp mp 
c/c, = Ki =p + ne 
(p+n)Q m) 
mp 
.). Ki = n + pe 
(p+n)Q m) (p+ n)q 
- ; t 
mp mV 
c/co = Ln + pe e 
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Figure 6—Partial mixing with piston flow and _ short- 

circuit: liquid from the zone of perfect mixing displaces, 

and is replaced by, an equal volume from the stagnant 
region. 


Plotting c/c, vs qt/V on semi-log paper, for the two 
regions of case A, gives curves such as shown in Figure 4 

(B) Part of the feed goes to the stagnant zone and 
displaces an equal volume by piston flow around the zone 
of perfect mixing, directly to the outlet, as indicated in 
Figure 5. The outgoing material is made up in part by 
this stream, in part by liquid forced out of the zone of 
perfect mixing and also by the part of the feed going in 
short-circuit. 


Again, two periods have to be considered, before and 
(1—m)V 


Pq 


after time: t = 


(1-—m)V 
< — 
Pq 


a)O<t 


nqcrdt = nqe’dt + mVdc’ 


de’ ‘ nq 


dt mV 


c’ =O 


But, 
c = nc’ + pc, 
_ nq 
mV 
.). C/Co = p tne 
(1 — m)V 
b) t 
Pq 
nq 
mV 
c’/Cs = € 
Now, 
c = nc’ 
— | 
mV 
~). C/Cy = ne 


Curves of c/c, vs a on semi-log paper, for these 


cases, are as shown in Fig. 


(C) A certain amount of liquid, corresponding to a 
fraction s of that part of the feed which enters the zone 
of perfect mixing, passes out of the latter into the stagnant 
zone. As a result an equal volume is displaced, in piston 
flow, from the stagnant section into the zone of perfect 
mixing, as shown in Figure 6. 
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[he two periods to be considered are before and 


: (1—m)V 
after time: t= - 
nsq 
(1—m)V 
ayo<te<d 
nsq 


Throughout this time period, liquid of composition 
c, is displaced into the zone of perfect mixing from the 
stagnant one. 


nqcrdt + nsqe,dt = nsqe’dt + nqe’dt + mVdc’ 


de’ nq emt nsqCo 
dt mV mV 
nq 
“ ~ mV oe ee 
‘sy " r 
c’/Co = + Kee 
°  s+1 
’ ° ’ . 
Since c’ = c, at time t = O 
5 1 
Ke = 
s+ 1 
Also, c = nc’ 
) ; n + — “(s+ 1) t 
-). C/o = ste 
/ ©o s+! 
(1—m)V 
Bt > 
nsq 


Liquid entering the zone of perfect mixing now has 
the composition c’ of the liquid displaced from the zone 
of perfect mixing into the stagnant one, at time: 








Figure 7—Partial mixing and short-circuit in a chemical 
reactor. 


Now, taking account of the reactant in the feed and 
exit streams as weil as that converted in the reactor, one 
can write (2°); 


PU iisn cava eee (5) 


nqcr — nqce’ = 


where: r’ = reaction rate, lbs./(ft.3) (hr.), in the zone of 
perfect mixing. 


mV = actual volume of liquid undergoing reaction. 
Equation (5) can be rewritten: 
; r’mV 
c’ = cr — 
nq 


Introducing this value of c’ in Equation (4) gives: 


(1— m)V oe 
t - — or, att — ¢). : r’m\ 
nsq CR “in iticscentanenvewas (6) 
q 
nqcrdt + nsqe’ ¢—;~,, dt = nsqe’dt + nqe’dt + mVde’ 
nq a m)V 
dc’ nq Consq mee es nsq ) 
RL - +e 
dt m\ (s + 1)m\ 
(s + 1) (1 — m) nq nq 
ty s? nsqt sm mV res ie mee me 
ise es ; © S3e 
(6. t)* ~ (e -F Ba 
(s + 1) (1 — m) 
‘ 1 s 1—m ain 
K; = - _— - le 
i ae i ie) m(s + 1) 
(1 +s) (1 m)~ _ nq 
sm mv (1 + aia: 


: ns? 4 n 1+ s 1 —m 4 nsqt 
c/o = a Pe = 
wae oe ae er i+s m mv } ‘ 


Curves for the last two cases are shown in Figure 6. 


If mixing was perfect in the entire reactor, the con- 
centration of the efluent would be cy and related to cy 


5. Design procedure for tank flow reactors when : ; 

partial mixing occurs by the following: 

As pointed out by Smith ), design procedures for aise de r\ eee F (7) 
tank flow reactors are still unavailable when partial mixing q 


occurs. The problem may now be solved, in general, by 
introducing the concepts already mentioned. 

For the particular case of partial mixing and short- 
circuit, as shown in Figure 7, under steady-state flow 
conditions, the incoming solution has a composition cy, 


where: r = theoretical reaction rate, Ibs./(ft.3)(hr.), 
based on concentration c),. 
Combining Equations (6) and (7) gives: 


r’ 
) 


c = cr — m(cr — Cy 
that in the zone of perfect ni... ing is c’, and the outgoing t 
liquid has a composition c. e e ) . 
e m = ; (8) 
A material balance at the outlet givac Cy eu / r 


(1 — n) qcr + nqe’ = qe 


or,c = (1 
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While Equation (8) is a general solution for partial 
mixing and short-circuit, many particular solutions can 
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be obtained depending on the order of the reaction taking 
place. 

Of all the possible solutions, only the following two 
are considered here: 

a) Zero-order reaction with partial mixing and short- 
circuit. 

For such a reaction, the rate is independent of con- 
centration, so that r = r’. 


Then, Equation (8) becomes: 
C} =< 
m= 
Cr — Cm 
Thus, one could either obtain the value of c for a 
given value of m or determine the effective volume by 
measuring feed and outlet concentrations and calculating 
the theoretical outlet composition. It is to be noted that 
in this case the short-circuit effect cancéls out. 
b) First-order reaction, with partial mixing and no 
short-circuit. 
In a first-order reaction the rate is proportional to the 
concentration of the reactant. If the concentrations were 


Cy and c’, the reaction rates would be respectively: 
r= key 


1 
F = KC 


Substituting in Equation (8), one obtains: 


Cy Cc CM 
mn = ; 
Cr — (msc 


Now, introducing the value of c’ given in Equation 


(4), 
noe (2 ad )( cm ) ree (9) 
Cr =- Cy ner —~— Ce + ct 


Since the case considered here is that where no short- 


circuit is present, n = 1 and Equation (9) becomes: 
fe = CM 
m = _— 
c Cr — Cm 
A 
r 
| 
5} 










ne 
. 


f 
oy 
a 
o en 
of Pia 
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Figure 8—Relation between actual and theoretical outlet 


concentrations for tank flow reactors when partial mixing 
occurs. 


110 


SAMPLE 


The Canadian Journal of Chemical Engineering, June, 1959 


which can be re-written: 


—— m (*) + (1 — m) 
Cc CM 


Cr . ° 
The value of “" can be obtained easily from the above 
. ; 


° ° ae . ° Cr 
expression, or from Figure 8, for given values of —* and m. 
£ £ Che 


The designer is thus in a position to predict the per- 
formance of a tank-flow reactor. 

For cases where short-circuit is present, where higher 
order reactions are involved or where there is piston 
flow, studies are actually under way which will be the 
object of a future publication. 


Apparatus 


Of the different cases studied theoretically, that in- 
volving partial mixing and short-circuit, as described in 
case No. 2, seemed to be of more immediate interest 
from an experimental point of view. An apparatus was 
set up, as shown in Figure 9, in which could be studied 
simultaneously the variations of effective volume and 
short-circuit with the level of agitation. 

The tank used had an internal diameter of 30”. A liquid 
depth of 30” was assured by a properly located outlet of 
}” pipe. The agitator was an experimental model by 
Chemineer, a 23” simple turbine with inclined blades 
being used. A flow regulator and a rotameter were 
inserted in the ?” feed line. The feed inlet was 12” from 
the surface of the liquid and the agitator itself was placed 
in between, at a distance of 6” from the surface. 


The agitator speed could be adjusted at the proper 
value through a variable speed reducer. 


Experimental procedure 


At the beginning of a run the reservoir was filled to 
the appropriate level with a solution of NaCl, approxi- 
mately 1/20 N. After running the agitator for some time 
at the desired speed to allow for the mixing pattern to 
be fully developed, water was introduced suddenly at the 
required rate of approximately 1.15 gallons/min. 

Samples were taken at regular time intervals of five 
minutes and the NaCl concentration determined volu- 
metrically with AgNO,, using K,CrO, as indicator. 


Results 


The experimental results obtained were plotted, as 
suggested already, on semi-log paper. Figure 10 shows the 
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Figure 9—Experir ental apparatus. 



























Figure 11—Effective volume of mixing vs. agitator speed. 
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Figure 10—Experimental results: relative concentration vs. 
time. 
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Figure 12—Short-circuit effect vs. agitator speed. Figure 13—Variation of m and n with agitator speed. 
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straight lines drawn through the experimental points at 
different levels of agitation, for values of agitator r.p.m., 
designated as N, ranging from 0 to 210. Values of m and 
n were determined from them for each agitator speed by 
the method described. 

Even at an r.p.m. of 0 given values of m and n were 
obtained; they were designated as m, and n,. For any 
given value of N, then, the value of m — m, might 
represent the contribution of the agitator r.p.m. to the 
level of mixing. 

It has been found that when plotting (m — m,) vs N 
on log-log paper, as shown in Figure 11, a straight line 
of slope 2 is obtained which can be represented by the 
equation: 

m — m, = a N? 


For the particular system studied, the constants m, 
and a were found to be 0.38 and 1.25 x 10° respectiv ely. 

The value of 0.38 for m, indicates the level of mixing 
obtained when the agitator is at rest, through the agita- 
tion produced by the introduction of the feed. Such an 
effect anyones to have been mentioned first by MacDonald 
and Piret *), when measuring mixing times in a vessel. 


As the value of m reaches 1, when the agitator speed 
is increased, a maximum value of (m — m,) is reached. 
This value remains constant thereafter even if the agitator 
speed is increased because the system behaves as one 
giving perfect mixing. 

Proceeding in a similar way, the values of (n — n,) 
vs N were plotted on log-log paper as shown in Figure 

2, giving a straight line of slope 3. The values of n are 

thus related to the agitator r.p.m. by the equation: 
an, = b N® 

where, for the particular system studied, n, = 0.77 and 

b= 2x 104. 


When plotting, on a linear scale, the values of m and 
n as a function of N, as in Figure 13, both variables are 
found to reach a maximum value of 1 at an r.p.m. of 
approximately 220, indicating that the system behaves 
then as one giving perfect mixing. 


Conclusions 

The method described to determine mixing efficiencies 
in continuous flow systems should prove to be a very 
useful tool in ev raluating the actual performance of 
mixers. Not only is it possible to obtain optimum results 
with a given mixer but the relative performance of mixers 
of different design can be established since efficiency 
measurements are possible. 

The contribution to the design of tank flow reactors 
should also be of importance because a method is now 
available to take into account the effects of partial mixing. 
Improved designs, with the best possible use of agitators, 
should result. 

The experimental results obtained to date are very 
promising. While they do not draw a complete picture 
as yet, due to the limited extent of experimental work 
which could be performed, they point to a number of 
investigations which should lead to a comprehensive 
knowledge of mixing. 


A new experimental set-up providing for a wider 
variety of investigations is nearing completion. Among 
others, the factors to be studied include the position, 
speed of rotation, shape and size of the agitator, the 
physical properties and the relative depth of the liquid, 
the location of the feed, the manner in which it is intro- 
duced, its rate and inlet velocity, and the outlet position 
relative to the inlet. 


Other types of agitating systems will be considered, 
including those with internal recirculation and of the 
draft-tube type. 


The study of systems involving chemical reactions is 
already under way. 


Nomenclature 


a = Constant, relating effective volume to agitator speed. 
b = Constant, relating short-circuit effect to agitator speed. 
c = Volumetric concentration, lbs./ft.* 
cr : concentration of the feed. 
Cm : concentration of the effluent if mixing was perfect, 
when a chemical reaction is present. 
Co : initial concentration. 
c’ : concentration in the zone of perfect mixing 
*, Ki, Kj, Ke, Ks : Constants of integration. 
= Fraction of the total volume which is perfectly mixed 
m, : value of m when the agitator is at rest. 
= Fraction of the feed entering the zone of perf »ct mixing. 
No : Value of n when the agitator is at rest. 
= Agitator speed, R.P.M. 
Fraction of the feed going to the stagnant zone. 
= Rate of flow, ft.3/hr. 
= Reaction rate, lbs./(ft.*) (hr.) 
Fraction of part nq of the feed entering the zone of perfect 
mixing which subsequently passes out into the stagnant 
zone. 
= Time, hrs. 
(1 — m)V 


nsq 


V Volume, ft.* 
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Sodium Hydroxide-Sodium Sulphate 
Stripping of Uranium from 
Amine Extracts 


A. I. BELLINGHAM? and R. SIMARD® 


Batch shake-out tests on two tertiary amines and 
one secondary amine to determine the effects of 
sodium sulphate concentration and pH on the distri- 
bution of uranium between 0.1N amine in kerosene 
and aqueous sodium sulphate solution showed that 
transfer of uranium from the solvent to the aqueous 
phase was favored by increasing pH values up to 5.5 
and by increasing the aqueous sodium sulphate con- 
centration up to 20% w/v Na,SO,,. 

Continuous counter-current stripping of tri-iso- 
octylamine (TIOA) extracts with 20% w/v Na,SO, 
was also carried out in a three-stage mixer-settler 
arrangement in which the pH was controlled by the 
addition of 10% NaOH to the three stages. Further 
addition of caustic to the separated aqueous strip 
solution to raise the pH to 7.0-7.2 yielded a high 
grade uranium precipitate (88% U,O,). The contin- 
uous testing confirmed the results of the shake-out 
tests and showed that the procedure gave results com- 
paring favorably with other stripping methods accom- 
panied by a lower overall reagent cost. 


A solvent extraction of uranium from sulphuric 
acid solutions produced in the leaching of ores has 
been the subject of much investigation in the last few 
vears. At the U.S.A.E.C., Oak Ridge National Laboratory, 
in particular, K. B. Brown, D. J. Crouse and others have 
made extended studies “') of long chain amines which 
have proven to be efficient extractants for uranium. 


The amine solvent extraction is best explained by an 
ion exchange mechanism which may be represented as 
follows: 


(RsNH)2 SO, + [UO.(SOx,) 2|= 


z:2=2 (R3NH)2 . UO(SO,) » + SO= 


It is well known that the loaded amine extracts can be 
stripped by a similar exchange using acidified solutions 
of chloride or nitrate salts, or w ith sodium carbonate, 
ammonium carbonate and alkalis which react irreversibly 
as neutralizing agents. 

In the course of studies being made on the application 
of solvent extraction to the hydrometallurgy of Canadian 
uranium ores, work relating to the use of certain amines 
has shown that the loaded extract can also be stripped 
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effectively by sodium sulphate solutions at a controlled pH 
favouring uranium in the aqueous phase. This led to an 
investigation of continuous counter-current stripping 
which has been incorporated in a pilot plant flowsheet. 

In this work, three amines were investigated by small 
scale batch testing. These were (1) tertiary tri-iso-octyl 
amine supplied by Carbide Chemical Company , Montreal, 
Que., (2) a tertiary amine labelled C,-C,, mixture sup- 
plied by General Mills, Inc., Kankakee, Ill., (3) a secon- 
dary dodeceny] amine labelled 9D-178 supplied by Rohm 
and Haas Company of Canada Limited, West Hill, Ont. 
In the batch testing all three products responded favorably 
to stripping with sodium sulphate solution. However, since 
the tri-iso-octyl amine, at the time the work was done, 
appeared to have some advantages over the other two in 
the extraction section of the process, it was chosen for 
the continuous testing. 


Procedure 

Batch cross-current stripping was done in open 
beakers with mechanical stirrers using 250 ml. of 0.1N 
amine (in kerosene) which had been previously saturated 
by contact with uranyl sulphate solutions at pH 1.4-1.5 
Three or four 50-ml. fractions of sodium sulphate solusien 
were contacted successively with the extract. For each 
stage of stripping 10% NaOH was used to obtain the 
required pH at equilibrium. The strip fractions were 
separated after each stage in a separatory funnel and 
assayed for U,O,. The ‘final stripped solvent was also 
assayed for U. O.. 


Two main variables to be studied were Na,SO, con- 
centration and pH. The uranium distribution curves, Co 
vs. Ca expressed as U,O, g/l in organic and aqueous, 
were plotted as a function of pH and sulphate concen- 
trations. From these data it was possible to predict the 
feasibility of an efficient continuous stripping circuit 
using three counter-current stages at a controlled pH. 
This was incorporated in a bench-scale extraction unit 
already treating uranium leach solutions. 


Results 
A. Batch cross-current stripping 

E oe data are represented graphically in Figures 
1, 2 and 3. The pH as shown in Figure 1 was found to 
be a csidieal variable. For efficient stripping or for less 
than 0.1 g U,O,/1 in the solvent, the pH should be 
controlled at 4.5 to 5.0. Higher pH values will cause local 


precipitation in the first stage where uranium concentra- 
tion is high (15-20 g U,O,/1). 
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Figure 1—Distribution curves. O.IN tri-iso-octyl amine vs 


20% w/v Na.SO, Effect of pH. 


The effect of varying the sodium sulphate concentra- 
tion is illustrated in Figues 2 and 4 for 0.1N TIOA (Tri- 
iso-octylamine ). In Figure 2 the dotted lines would repre- 
sent the true equilibrium if no sulphate were picked up 
in the first cross-current stage. At this stage, for a volume 
ratio of solvent to strip of 5 to 1, the Na, SO, increase 
in the strip is approximately 3%, so that the curves for 
water, 10% and 20% of Na,SO,, are in fact, 3%, 13% 
and 23% Na,SO, respectively. This explains the apparent 
dip in the curves of Figure 2. The data from the corrected 
curves of Figure 2 were used to plot the distribution 
coefficients of Figure 4. 

In view of the results obtained with tri-iso-octylamine, 
other amine extractants were studied. Equilibrium data 
for Rohm and Haas secondary amine 9D-178 and for 
General Mills tertiary amine C,-C,, are represented in 
Figure 3. 


B. Continuous counter-current tests 


To confirm the equilibrium data obtained from batch 
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Figure 3 — Distribution curves. pH 4.75-20% ‘w/v Na:SO,. 
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strip at pH Effect of Na.SO, concentration. 


testing, a series of continuous extraction and stripping 
runs was carried out in a small pilot plant extraction unit 
designed to test various leach solutions obtained from the 
treatment of Canadian uranium ores. 

The flow diagram is shown in Figure 5. Figure 6 is a 
detailed sketch of the 4-stage mixer-settler unit currently 
used for extraction. The three-stage stripping unit is a 
scale-down of the extractor. The two units operating in 
a closed organic circuit were designed for a throughput 
of 400-600 ml. feed solution per minute. Loaded extract 
from the extractor flows counter-currently to a stream 
of Na,SO, solution through three mixer-settlers. The pH 
of the aqueous phase in the first mixer is controlled from 
a pH recorder-controller which actuates the sodium 
hydroxide valve. Stage 2 and 3 are controlled manually 
to maintain the set pH as read by a Beckman pH meter. 

The operating conditions and results of a number of 
continuous runs are shown in Table 1. Data from run 
532 were used to establish the operating equilibria of 
Figure 7. 






Organic Solvent = O1N TIOA 
pH of Strip = 5 0 


50 


45 
40 
35 
Ca 
— 30 
Co 


° 2 4 6 8 10 (2 14 16 18 20 
% W/V NazSQq IN STRIPPING SOLUTION 





Effect of Na.SO, concentration on distribution 
coefficient. 
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TABLE 1 


ConTINUOUS EXTRACTION AND STRIPPING TESTS 
OPERATING CONDITIONS AND RESULTS 





TABLE 2 


COMPARATIVE RESULTS OF VARIOUS STRIPPING METHODS 


| Reagent Precipitate Recycle 
Run No. 529 530 Strip Solution Consumption Grade Barren Strip 
— -—— - lb./Ib. U308 % U3;0s Treatment 
Flows ml/min —— ~—-- 
Pregnant Feed Soluticn a 305 Sodium Sulphate 
Organic: 0.1N TIOA, 100 20% Na2SO, 
Strip Solution: NaOH precipits ition | NaOH-1.0 85-88 None 
20% W/ ly NasSO, ‘ 38 ae - _ 
Silinn Chloride 
Ass says 1.0M NaCl NaCl-3.0 Removal of 
Pr gn int Solution, pH 1.60 | 1.6 1.60 0.05M H.SO, H.SO,-0.25 excess sodium 
U;0s gm, fH fas | 4a 1.81 NaOH precipitation | NaOQH-—0.70 sulphate by 
Raffinate Stage 1 a 1.75 1.7 0.98 crystallization 
‘ i sia 1.85 £5 0.028 — 
7) <3 = i 0.0: ‘ 0.001 Sodium Carbonate 
Se = 0.075 0.002 | <0. <0.001 10% NasCO; NaOH-1.30 Removal of 
Organic Extract oA9 Ko | 5.46 NaOH precipitation excess sodium 
Strip Stage 1 | 13.75 | 14.5 | 12.30 sulphate by 
" ie 5.39 | 3.3! 2.92 crystallization 
v = oe 1.98 | 0.6 0.41 and conver- 
Recycle Solvent 0.48 | 0.03 sion of excess 
Oo 30 | NaOH to 
: Average operating pH of | carbonate 
ine Strip Stage 1 | 4.4 | using a Car- 
Be ie | 4.2 | bonation 
. ae | 42 | tower. 
NaOH Consumed | | 
ng For strip control and precipi- 
: A, poe | 
“a tation Ib NaOH/Ib U;Os 1.1 by adding NaOH to pH 7.0-7.2 produced a high-grade 
Tne ; ; - 2 2 
Precipitate Assay U:Os% 80.5 product that could be settled readily and filtered well. 
The filtering rate on the settled portion was estimated at 
sa > 5 lb. dry precipitate per sq. ft. per hour. . 
™ ; : Reagent consumption for a and precipitation 
From the above table and from Figure 7, it can be seen 5 E f 
7“ BS wauee. § nary ’ ‘ : in all tests was close to 1 lb. NaOH per lb. U,O, and 
he that theoretical equilibrium conditions can be approached this confirms the theoretical consumption based on the 
in continuous operation. The pH control did not prove U,O,/SO, ratio in the extract. From the assays of the 
ut : a 
: as difficult as anticipated. After some preliminary exper extract the mol ratio in the solvent is close to 3 SO, to 
oe ience, 6 Secor oo. — be — to — 1 U. On this basis, complete neutralization to pH 7.0 
at a pri OF 4./7-U.1) which Is satisiactory 10F an On-o would be as represented in the following reaction: 
yH controller, At times there was some slight precipitation 
ym of uranium due to localized high pH at the point of 2 (Rs NH), . UO(SOs)s + 14 NaOH — 
im addition, but it redissolved and no trouble was exper- ; mes 8 Rs - a Na2U:0; + 6 Na:SO, + 11 H,0 
ly ienced. The phase separation in the strippers was rapid or 21 mols NaOH per mol Us ; 
. or close to 21 x 40/842 = approximately 1 Ib. NaOH /Ib.U;0, 
er. with the primary break occurring in 50 seconds. In the 
of mixing step of the stripping operation, the minimum A simple explanation for the mechanism of sulphate 
ai retention time so far used has been 5 minutes per stage stripping at controlled pH has been to assume that the 
ind has been adequate. uranyl ion forms a double salt with Na,SO, thus ogee. 
of ads eran, ; ; 
The precipitation of uranium from the sulphate strip sing ‘the activity of the anion complex [U O,(SO,),}4 
NaOH NazS0q Pregnant 
Solution Strip Storage 
Solution 
i Flowmeter 0 Flowmeter 
| 
pH Electrodes 
K | 
\ | 
Manual | Recycle Solvent 
EXTRACTORS 
Organic 
Hold- up 
Metering 
Pump 
High Grade Strip 
on Figure 5—Flow diagram. Extraction and sulphate stripping 
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Figure 7—Graphical representation of 
continuous counter-current extraction 
and stripping. 
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Conclusion 


The results of the test work reported here were suffi- 
ciently encouraging for the process to be tested on several 
uranium bearing solutions from the leaching of various 
ores. This additional work confirmed the work described 
in this paper ‘*-% 4) and demonstrated that sodium sul- 
phate stripping of uranium from amine extractants was 
simpler and cheaper than by chloride or carbonate solu- 
tions. The method is simpler in that the barren. strip 
solution may be recycled directly without any form of 
regeneration or purification whgreas other strip solutions 
must be retreated before recycling. The sulphate strip 
process, however, would necessitate more stringent pH 
control and instrumentation than other stripping methods. 
Table 2 shows a comparison of three common methods 


for stripping uranium from amine extracts. 
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Figure 6—Four stage mixer settler 
(14” lucite). 
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Nomenclature 

Ca = concentration U;Os g/1 in aqueous solution 
Co = concentration U;Os g/1 in organic solution 
Ca 

Co 

A = aqueous feed solution flowrate, ml/min 


= distribution coefficient 


O = organic solvent 
S aqueous stripping solution flowrate, ml/min 
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A Note on Thermodynamic Consistency 





of Ternary Vapor-Liquid 


Graphical and numerical integration methods are 
proposed for testing the thermodynamic consistency 
of ternary vapor-liquid equilibrium data, in which the 
integration of the Gibbs-Duhem equation can be per- 
formed under isothermal and isobaric conditions. No 
data-fitting procedure is involved in the proposed 
methods, therefore the uncertainty as to whether the 
deviations are due to experimental errors or to the 
inadequacy of the employed equation can be avoided. 
When one part of the data is less precise than the 
others, the proposed method may be employed to 
detect the inconsistent part. 


HERMODYNAMIC consistency has been the most reliable 
"hae of vapor-liquid equilibrium data. Especially for 
the ternary systems, the experimental difficulties involved 
are such that a preliminary set of data has to be tested 
against thermodynamic consistency before one can even 
trust his experimental technique. In the literature, many 
test methods have been proposed. The experimental data 
are tested either by fitting them with some integrated 
forms of the Gibbs-Duhem equation, or by graphical 
integration of the Gibbs-Duhem equation along a par- 
ticular path. When data fitting procedure is involved, 
using equations such as those proposed by Benedict et 
al.) and Wohl), the uncertainty as to whether the 
deviations are due to experimental errors or to the inade- 
quacy of the employed equation can not be avoided. 
Methods which involve graphical integration along a 
particular path may also encounter practical difficulties. 
Recently Herrington (3) has “er two methods, one 
of which involves integration of the Gibbs-Duhem equa- 
tion along a line parallel to one of the bases of the 
triangular coordinates, thus keeping the mole fraction of 
one component constant. The other method, which modi- 
fies the method of Colburn and Schoenborn ‘*), involves 
integration of the Gibbs-Duhem equation along a straight 
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line connecting one of the vertices to a point in the 
opposite base. In this method, the ratio of two of the 
mole fractions is kept constant. He also proved that these 
methods are applicable to systems with more than one 
liquid phase, provided that the whole liquid is treated as 
a single phase by using the overall mole fractions and 
activity coefficients, etc. The integration can therefore 
actually go through the two liquid phase region without 
any discontinuity. The purpose of keeping one of the 
mole fractions or the ratio of two of the mole fractions 
constant is to simplify the integration of the Gibbs- 
Duhem equation. However, to restrict some relation 
among the mole fractions has several disadvantages: 


1. Since most of the vapor-liquid equilibrium data are 
measured at constant pressure, a restriction in the mole 
fraction will force the integration to go through a non- 
isothermal region, to w hich, strictly speaking, the Gibbs- 
Duhem equation derived from isothermal and isobaric 
conditions cannot be applied. 


2. Unless the experimental points lie exactly on the 
straight line required for these methods, some interpola- 
tion or extrapolation may be involved which will cause 
additional uncertainty. 


3. These methods may not go through the region of 
more reliable or more abundant data at ease without some- 
times covering some undesirable region. 

The purpose of this investigation is to introduce a 
method in which the integration of the Gibbs-Duhem 
equation can be performed along any curve in the tri- 
angular coordinates, especially along the curve of constant 
temperature and pressure w there the method is thermo- 
dynamically exact. 


Integration of the Gibbs-Duhem equation 

Under isothermal and isobaric conditions, the Gibbs- 
Duhem equation for a ternary liquid solution may be 
expressed as 

x1 d log ¥1 + x2 d log yz + x3 d log y: = O (at T,p)...(1) 
where x,, x, and x, are the liquid mole fractions of 
components 1, 2 and 3 respectively, and y,, y, and 


are the liquid activity coefficients of components 1, 2 
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and 3 respectively, using pure liquids as standard states. 
The Scatchard’s excess free energy of mixing may be 
expressed in the form ‘°) 
Q =) gi x2 log Y2 + x3 log Ys 
Then, it is seen 
dQ = logy: d xi + log y2d x2 + log y3d x3 (at T,p). .(3) 
Since Q = 0 at the three vertices of the triangular co- 
ordinates, an integration of Equation (3) along a curve 
connecting two of the vertices will give 
c =0 x.=0 x.=0 
- 5 


e 


rc =0 x; =0 x.=0 


log yi d xi + | log y2dxo + 


Equation (4+) may be applied to ternary vapor-liquid 
equilibrium data obtained under isobaric conditions, pro- 
vided that the change of temperature along the curve is 
not large. For exact correlations, the integration can be 
performed along the curve of constant boiling temper- 
atures, which gives 


b b 
e - 


Qv — Qa = | logyidxi + | logy2dx2 + | log y3dxz3 (at T,p).(5) 


~ ~ 2 
a a a 


where a and b are any two points on the same curve. 
The problem is now how to evaluate these integrals. A 
direct method (Method 1) would be integrating each of 
the integrals graphically. This can be done by plotting 
log y, against x,, log y, against x, and log y, against x,. 
The points have to be taken from the curve in the 
triangular coordinates in the sequence from point a to 
point b. The same sequence has to be used in each plot. 
Since 


dx; + dx» + dx 


b 


Qpn — Qa = ioe 1 - ae : i oe OOD 
J Y3 Y3 


a a 
Therefore, only two plots, namely log (y,/y,) against x, 
and log (y./y,,) against x, will be necessary. If a and b 
coincide to one point, that is if the integration is taken 
along a closed curve, or a and b are two of the vertices, 
then the value of the sum of the integrals is zero, A some- 
what similar approach has been proposed by Broughton 
and Brearley ‘*) for testing isobaric data. However, instead 
of testing data along constant temperature curves, they 
assumed that the product t log y is independent of 
temperature. The validity of this ‘assumption seems ques- 
tionable even for binary systems 7). 

In some instances, calculations using relative volatilities 
are more advantageous than using the activity coefficients. 
If the integration is along the curve of constant boiling 
temperature and if the vapor phase can be assumed ideal, 
Equation (5) becomes 


V1 y2 
Ax: log + Xz log 


X) X2 


where a,., a,, are the relative volatilities defined by 





and A is the difference of the quantity in the parentheses 
between the points a and b. Should the curve be a straight 
line, the ratio of dx,: dx,: dx, is a constant, only one 
integration is needed. In this case the right-hand side of 
F.quation (6) becomes, 


b 


( Xo — Xa 2 
| | log a oa a — . log y dx;. 
vs Xib — Xia ¥3 


The evaluation of the integrals of Equation (5) may 
also be carried out by means “of a numerical integration 
(Method 2). This method has the advantage that it can 
a applied to almost any set of data without using any 
graphical procedure, provided that the experimental 
_ - reasonably spaced. Let the points be a, 1, 2, 

. 9, b. Assuming that all the curves of log y against x 
can approximated by segments of straight lines con- 
necting the points, the first ‘integral of Equation (5) can 
be calculated by 


1 
log y1 dx; = ; (log Yia + log Yi1) (Xin. — Xta) + 


(log yn + log yi2). (X12 — xu) +... 


a 
Xip log Yi - Yi9 — Xia log Via Yu + as 
da 
Xi log (Yia/Vi2z) + X12 log (Yu/Yis) #2. + ve 
dy 
dc 
X19 log (Yis/¥i») (10) : 
If 
a | 
Ww wae X;, 1s the mole fraction of component 1 at point a, in 
Via is the activity coefficient of component 1 at point a, 
etc. Similar expressions may be obtained for the second w 
and the third integrals. Combining these expressions to- W: 
gether with the Quantities Q, and Q,, which are defined co 
by Equation (2), the followi ing equation is obtained: ze 
Via Via il 
Xia log 2 Xi log + Xi log : + + 
Yu Yi Yi3 
X19 log 
Xyy log . + xy» log 
Y 3b 
This summation is simply that for each component, 
multiply the mole fraction at any point with the logarithm 
of the ratio . the activity coefficients of a point preced- 
ing it and of a point following it, except the first and 
last points, w ae the activity coefficient of its own is . 
used if there is no point preceding it or if there is no 
point following it. It is convenient to tabulate the calcu- Fi 
lation as shown in Table 1. eq 
If a and b coincide to one point, that is, if the curve 
is closed, then every point has a point preceding it and 
T 
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TABLE 


ACETONE-METHANOL-WATER AT 100°C. (4) 


(log ¥1p — 
log Vir) 


—0.01406 
.158 ri —(0.01506 
.119 .5035 —0.01297 
.089 bh —(0.01604 
.066 .623 —0).00283 
.028 3! ). 00088 
.033 5 -(). 00378 
043 a 00336 
.069 con -0.00128 


0,100 4925 00211 


= = (0.12284 
> — = 0.18122 
Repeat the above calculation without run Nos. 29, 
= + = 0.12051 
> — = 0.12767 


a point following it. The summation is to be performed 
as if every point is a point inside a curve. Experimental 
data over the entire concentration range may be con- 
veniently divided into several small sets and the thermo- 
dynamic consistency can be tested on each of them. In 
doing so, the less reliable region may be easily spotted. 
If one of the experimental points is much in error in 
a set of data, it may be singled out by successively includ- 
ing the experimental point, one by one, in a set. 

An example ef this calculation is shown in Table 1, 
where the data for the ternary system, acetone-methanol- 
water at 100°C. ‘*) are tested. The values of the activity 
coefficients are taken directly from the original article ‘). 
These activity coefficients are the thermodynamic activity 


Figure 1 — Thermodynamic consistency of vapor-liquid 
equilibrium data for the system acetone-methanol-water at 
100°C.(@). Percentages indicated are overall deviations. 
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log Y2 x» (log Y2, 


0.0398 


1 


= log ¥3 x; (log Y3p — 
log Y21) log Y2:) 


+0.00316 : 0792 —0.00136 


0560 -(). 00200 ‘ 0810 +0.00054 
0170 00536 762 0785 —( 

0110 00665 7 0828 +0 

0678 01747 5 0588 

1358 -0. 00888 0382 01533 
1166 00590 763 0394 02258 
1069 01140 ‘ 0678 02938 
0660 01486 oe 0795 02082 


0318 -(). 00659 7 0962 01179 


Total = — 0.05838 
©) Deviation = 38.4% 


Total = — 0.00716 
©) Deviation = 5.8%; 


coefficients for the liquid, and are not based on assuming 
an ideal gaseous solution in vapor phase. Since the pressure 
effect on liquid activity coefficients is usually small, it is 
neglected in this cance For the ten experimental 
points (Nos. 29, 31, 33, 35, 38, 39, 37, 36, 34 and 32) which 
form a curve as he fai on the triangular diagram 
(Figure 1), the overalf deviation in thermodynamic con- 
sistency is aaa. 38¢-, while the experimental points in 
the middle of the triangular diagram (Nos. 14, 17, 16, 18, 
21, 23, 26, 28, 30, 25, 24, 20 and 19) indicate a deviation 
of 10°-. The deviation is defined as 


se 
(ae “ ) Cae 
~>+)-(2 


( 
) 


j 
Aa 


f \ 
204— —4- —)80 


Figure 2— Thermodynamic consistency of vapor-liquid 

equilibrium data for the system n-octane-ethylbenzene- 

cellosolve at 760 mm. Hg. pressure (’). Percentages indicated 
are overall deviations. 
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It is obvious that the consistency of the data in the middle 
of the triangular diagram is better than the consistency of 
the points Nos. 29 . . . 32. If the calculation as shown in 
Table 1 is repeated, excluding the experimental points 
Nos. 29, 31 and 33, the ©, deviation reduces to 6%. It is 
seen that these points seem to be inconsistent. 


The application of the numerical integration method 
to a ternary system at constant pressure is shown in 
Figure 2, in w hich the data of the system n-octane-ethyl 
benzene—cellosolve ‘*) are tested. The experimental points 
are grouped in such a way that the boiling point differ- 
ences are kept at a minimum. It is seen that the deviations 
are not uniform throughout the whole concentration 
range. The experimental data on high cellosolve concen- 
trations seem to be less consistent than the data of low 


cellosolve concentrations. 


Although this numerical integration method is suffi- 
ciently accurate for experimental data of average precision 
in most cases, it can only be reliable if the points are not 
too few and are close enough to warrant the segmental 
line approximation of the log y vs. X curves . For fewer 
and widely separated points, a better approximation which 
utilizes the parabola approximation of any set of three 
consecutive points may be used (Method 3). This is 
similar to the Simpson’s method of aad integration. 
For data of very high precision, this method is also 
needed to detect any less precise points. 

For any set of three consecutive points, say x,,, Xj. 
and x,. (mole fractions of component 1 at points 1, 2 and 
3 respectively ) with values of log y,,> log y,, and log 
Y.3 uf 8 parabola is drawn to pass them, the area under 
the curve of log y against x is equal to 


3 
( 1 
log 1 dx; = 3 = ! ) log y1 it = (x13 — X11)”. 
) 6 


log (Yi2 Vi1) 
C12) 
X} 1 ae = Xa 


Similar expressions may be obtained for the second and 
third integrals of Equation (5). In Equation (12) the first 
term on the right hand side is siniles to that of E quation 
(10). Therefore, the second term serves as a correction 
to the second method and the order of magnitude com- 
pared to the overall summation of Equation (11) will 
indicate whether the third method is necessary or not. 
In a practical calculation, this method requires a set of 


odd number of experimental points. The sum of all the 
integrals should add up to (Q, Q,) for any set of 
points which constitute a curve in the triangular co- 
ordinates. 

Still more accurate numerical method does not seem 
to warrant discussion. Only some points of interest will 
be mentioned here. First, since it is not necessary to 
restrict any relation between the mole fractions, the inte- 
gration does not have to go through a multi-liquid phase 
region where there are almost no actual experimental data. 
The integration can stop at a point on the saturation 
curve or go through along a tie line. The integration 
along a tie line is not necessary since the terminal Q 

values can be easily calculated. Instead, it is recommended 
to integrate along the saturation curve where the experi- 
mental data are usually available. Second, it is customary 
to plot the vapor- liquid equilibria on a triangular coordi- 
nate paper using the vapor composition, y, of one of 
the components as parameter. In this case, the Gibbs- 
Duhem equation can be modified to become: 
Pr ye Pi Vz 
d log pi° = x2 d log — + x; d log Pat (13) 
pe Ps 
where the vapor phase is assumed ideal and the vapor 
composition of component | is a constant. Integration of 
Equation (13) along the line of constant y, will be 
further simplified if the temperature variation is small: 


b b 
r 


| xed log y: + | xsd logy: =0...... (14) 


a a 


Numerical integration can be performed using either the 
second or the third method by tabulating x, versus log y, 
and x, versus log y,. 
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Some Thermodynamic Properties 


of the System Formic Acid-Water'’ 


A. C. PLEWES?, 


The integral latent heats of vaporization of for- 
mic acid-water mixtures at 760 mm. Hg. are reported 
and used to prepare an enthalpy-concentration chart 
for the system. The data provided a basis for the 
estimation of the heats of mixing in the vapor phase, 
and the results indicate that the exothermic heat of 
mixing may be attributed to the partial dimerization 
of monomeric formic acid coupled with the associ- 
ation of water and acid molecules. 

Vapor-liquid equilibria data for the system are 
also reported. 


bers paper is one of a series“: *) reporting the results 
of studies carried out at Queen’s University on the 
properties of liquid mixtures in an effort to obtain useful 
data for industrial design projects and to explain the 
behavior of two unlike molecules in homogeneous solu- 
tion. The experimental work involved the measurement 
of the integral latent heats of vaporization of formic 
acid-water mixtures. Vapor-liquid equilibria data at 760 
mm.Hg. were also obtained as complete temperature- 
composition data were required in the preparation of an 
enthalpy-concentration chart. The previously reported 
equilibria data of Othmer ‘) could not be used as no 
temperatures were recorded and the data were obtained 
at an average total pressure of 750 mm. Hg. Similar studies 
of acid-water systems involving other homologous 
organic acids are being continued. 


Experimental equipment and procedure 

The experimental heats of vaporization were obtained 
in a modified version of a calorimeter used by Dana (. 
The calorimeter (Figure 1) is composed of a heater, C, 
and a silvered, evacuated Dewar flask, B. B is surrounded 
by a concentric glass cylinder, $, which is sealed to the 
top of the Dewar flask. Liquid from a constant head 
device enters S through Y and is boiled a controlled 
pressure of 760 mg.Hg. The chamber, S, is connected to 
a reflux condenser at R. Solution at the bubble point 
enters A through the capillary, D, in the base of the 
Dewar flask. The solution entering A is evaporated as 
quickly as it enters, thus ensuring evaporation at constant 
composition. The vapor produced in chamber A is with- 
drawn through the sidearm N, condensed, and retained 
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Figure 1—Calorimeter. Figure 2—Electrical circuit. 


for future sampling and analysis. The vessel, S, is heated 
externally by a chromel he: ating element imbedded in 
insulation and acts as a preheater for the liquid feed. 
Heat transfer between chamber A and the surroundings 
is reduced to a negligible amount by the evacuated 
Dewar flask and chamber S$. For the system formic acid- 
water the maximum temperature difference at steady state 
between chambers A and S$ was less than two centigrade 
degrees. The — ns consists of an electrical resistance 
element immersed in a bath of high- boiling silicone oil. 
The resistance aces and bath are completely enclosed 
in a glass tube. With this modification in the original 
calorimeter design corrosive liquids may be studied. 

A diagram of the electrical circuit is shown in Figure 
2. A 24-volt battery served as a power source and was 
connected in parallel with an 18-volt battery as ballast. 
This arrangement reduced the current fluctuations in 
the heater. The specific gravity of the battery electrolyte 
was kept within the range 1.20-1.30. The current was 
obtained from measurements of the potential drop across 
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Figure 3—Integral heats of vaporization of formic acid- 
water mixtures. 


a standard resistance of 0.5014 ohms. Two 10,000 ohm and 
one 1,000 ohm standard resistances, connected in parallel 
with the heating element, were used for voltage measure- 
ment. With this arrangement the power input to the 
heating element could be measured within = 0.2%. 

The vapor-liquid equilibria data were obtained in a 
still designed by Fillis ‘°). Ellis has tested this still for a 
number of binary mixtures involving both high and low 
activity coefficient systems composed of close and wide 
boiling components and compared the results with those 
obtained in stills of the Othmer, Fenske and Gillespie 
types. On the basis of activity coefficient tests he has 
shown that better experimental data are obtained with 
this new design. 


Analytical formic acid (90 wt.% acid) was obtained 
from the British Drug Houses and distilled in a column 
packed with | in. glass helices to a depth of five feet. The 
operation of the column at atmospheric pressure and a 
reflux ratio of 20:1 gave a middle third distillate of 
99.6 we. acid. This material was subsequently frozen 
and distilled under vacuum at room temperature to 
remove the residual water. The density of the final 
product was 1.2218 gm./ml. at 20°C. Once distilled city 
vater was used without further treatment. 


I he condensat¢ salniples were titrated Ww ith O.IN 
barium hydroxide protected from atmospheric carbon 
dioxide. Phenolphthalein was used as an indicator. 


The calorimeter was tested by measuring the heats of 


vaporization of formic acid and water at 760 mm.Hg. 
Ihe mean value of three determinations of the heat of 
vaporization of the purified formic acid was 118.0 cal./gm. 
4 value of 117.8 cal./gim. for 100¢, acid was obtained by 
( strapolation of the heats of vaporization for acid-water 


mixtures ranging from 92.5 to 99.6 weo% acid. A value 


{ 


of 120.0 cal./gm. obtained by Brown is given in the 


International Critical Table \n experimental value 


of 115.4 cal./gm. 1 reported by Coolidge Using the 


Clapeyron equation and hi experimental vapor pressure 
und vapol 
value of 114.0 cal./gin. A value of 115.6 cal./gm. is listed 
500 of the National Bureau of Standards ‘). 


It is beheved that the data reported in this paper are more 


und liquid density data ¢ oolidge estimated a 


in Circula 


accurate than the ilues just listed. This opinion is 
pronipt ib the design ind operational details of the 
alorimeter and on the fact that the ext apolated value of 
117:8 i./gin. checked the results obtained with purified 
fori i Ihe « sperinental heat of vaporization ol 

irc) i 34.5 il./gin hich checks with the value ot 
539.4 cal./gm. given by Keenan and Keyes 
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Figure 4— Enthalpy-concentration chart for the system 
formic acid-water at 760 mm. Hg. 


[he isobaric integral heats of vaporization of formic 
acid-water mixtures were obtained over the entire range 
of concentrations. Two runs were made at each con- 
centration and an average value is reported. Runs were 
continued for a period of one to two hours to ensure 
steady state operation during the experimental measure- 
ments. The criterion for steady state was the constancy 
and similarity of feed and distillate compositions. The 
calorimeter was operated at a boil-up rate of 3-5 gm./min. 
in Chamber A. At this boil-up rate there was no entrain- 
ment and a determination could be made in a reasonable 
time. 

The vapor-liquid equilibria data were obtained by 
operation of the Ellis still for periods of two hours. A 
second series of thirty-minute runs was carried out in 
the 90-100 wt.% acid range when it was discovered that 
formic acid tended to decompose when mixtures of high 
wid content were subjected to prolonged refluxing. The 
decomposition Was found to be negligible if the 
refluxing time was less than one hour. Vhe still was tested 
with the system acetic acid-water to see if equilibrium 
conditions could be reached in thirty minutes. The data 
obtained checked well with published values ('. Ac- 
cording to Ellis) 20-30 minutes is sufhicient time to 
reach equilibrium. 


Experimental results and discussion 


Enthalpy-concentration data for the system formic 
acid-water are presented in Table | and Figures 3 and 4 
The integral heats of vaporization listed in Vable 1 were 
read from the smoothed curve through the experimental 
data of Figure 3. The specific heat and heat of solution 
data for formic acid-water mixtures reported by Glago- 
leva (|! ind the temiperature Composition data and 
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integral heats of v aporization obtained experimentally i in 
this study were used in the preparation of the enthalpy- 
concentration chart (Figure 4). 


PABLI 


EN THALPY-CONCENTRATION DATA FOR THE SYSTEM FORMIC 
\cip-WaATER AT 760 MM. HG. 


Integral Enthalpy of | Enthalpy of 

Heat of Saturated Saturated 
Vaporization, Liquid, Vapor, 
wt. % acid cal./gm. cal./gm cal./gim. 


Composition, 
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The enthalpy of each pure component, in the liquid 
state, was assumed to be zero at 0°C. The datum level for 
liquid formic acid was chosen as 4.226 cal./gm. at 8.3°C. 
and 760 mm.Hg. This corresponds to an enthalpy of zero 
for liquid formic acid at 0°C. and 760 mm.Hg. as the 
normal freezing point of formic acid is 8.3°C. Vapor- 
liquid equilibria data for the system are also included 
in the enthalpy-concentration chart. 

Liquid isotherms were estimated at 20 degree intervals 
using the heat of solution and specific heat data. The 
saturated liquid line was drawn so that it intersected the 
liquid isotherms at known temperatures and compositions. 
The saturated vapor line was located by adding the 
integral heat of vaporization at constant pressure to the 
enthalpies of the liquid line. Isotherms in the vapor 
region are not included in the enthalpy-concentration 
chart due to the lack of experimental heats of mixing 
in the vapor phase. 

In the final phase of the study the heats of mixing of 
formic acid and water in the vapor phase were estimated. 
The results provided information on the possible struc- 
ture of the vapor and reasons for the departure of the 
system from an ideal solution. The procedure used was 
similar to that described in detail in a previous paper 
for the svstem acetic-acid water ‘2). The method is 
reviewed briefly as follows: 

(1) The heats of mixing were estimated from the 
difference between the actual enthalpies and the — 
enthalpies for various mixtures. The ideal enthalpy of ; 
mixture was obtained by the addition of the enth: iain 
of the pure components at the ideal conditions of zero 
pressure and the temperature of the system. 

(2) The variation of enthalpy of the pure compon- 
ents with pressure was obtained by graphical integration 
of the thermodynamic relationship “ 


5V\ 
. dp 
ol a 


The ideal enthalpies of water vapor were 
by extrapolation of the data in the steam tables ‘” to zero 
e in enthalpy 


AH = H, 
p 


estimated 


pressure. In the case of formic acid the chang 
with pressure was determined down to 25 mm.Hg. by 


graphical integration of Equation (1) using the P-V-1] 
deta of Coolidge ‘'*). From 25 mm.g. to zero pressure 


it was assumed that the acid vapor was composed of 
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monomer and dimer, and that any resulting enthalpy 
change was due to association or dissociation reactions. 
Thus, 
AH = Ho;mm.H¢g. H,° = aAHeg 
where 
AHex = molal heat of the association reaction 
a = fraction of monomer dimerized at equilibrium 
; 1 : 
HCOOH ; , (HCOOH )» 


( Pdimer pe 


Pmonomer 


-tx,/ 1 acs 
4PK3} +1 


P = total pressure 


1 -4/ I 
4PK? +1 


as 1 


where 


~ 


| atte aa 


AH = aAHr 


since 


The change in enthalpy due to om dimerization 
reaction was evaluated using Equation (2) and the heats 
of dimerization and equilibrium constants Geena by 
Waring “'5). The total change in enthalpy with pressure 
for formic acid was then obtained from the sum of 
Equations (1) and (2) and the results are presented i 
Table 2. 

The estimated heats of mixing and dimerization are 
compared in Figure 5 and Table 2. The heat of mixing 
in the vapor phi ise is exothermic and approximately 906, 
of this heat effect may be attributed to the disnevbauatins 
reaction. At low acid concentrations the difference be- 
tween the heat of mixing and the heat effect due to 
dimerization is small and increases with acid concentra- 
tion to a maximum in the region of the azeotrope (77.4 
wt.G%acid ). This increase is believed to be due to the 
association of formic acid and water molecules. If it is 
assumed that the heat of mixing is due to (1) the heat 
of association of formic acid monomer to dimer and (2) 
the heat of association of formic acid and water molecules, 
the structures of the associated molecules in the vapor 
phase may be pictured as follows: 


Formic Acid Dimer 
i oO 


\ssocia ted Formic \cid Water Molecule 
O tH 


0 tl Oo tl 


Krom the above structures it is e\ ident that the heat effect 
pel associated molecule for the dimerization reaction is 
approximately twice that for the formic acid-water asso 
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TABLI 


HEATS O1 


Enthalpy ol 
Formic 
Ac id Vapor 
(zero pressure 
and dewpoint), 
cal. 


Enthalpy of 

Water Vapor 

(7erTo pressure 
and dewpoint ), 
cal 


Composition, 


( cal. 


wt. © acid em. em. 

- 0 642.9 642 
10 643.5 287.0 607 
20 644.1 288.0 572 
30 644.5 88.7 537 
40 644.9 289.5 502 
50 645.3 290.4 467 
60 645.7 291.3 433 
70 646.0 291.8 398 
80 645.9 291.8 362 
90 645.5 291.0 326 
100 286.1 286 

ciation. Unfortunately, there is no direct means of 


estimating the heat effect due to the association of formic 
acid and water molecules in the vapor phase. Other 
associated structures are possible, e.g. water with formic 
acid dimer, but the degree of association in more com- 
plicated structures than those considered should be quite 
small and the corresponding heat effects negligible. 
Vapor-liquid equilibria data required in the con- 
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Figure 5—Heats of mixing in the vapor phase. 
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MIXING OF ForMic ActpD-WATER MIXTURES IN THE VAPOR PHASE AT 760 MM. Ho. 


Enthalpy of 
Vapor Mixture 
(zero pressure 
and dewpoint), 
gm. 


9 
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Heat of Mixing D 
Enthalpy of ~ ee + 
Vapor Mixture a 
(760 mm. Hg. Hiatm, — Hpeo, Heat of a 
7 . . ae . oO 
and dewpoint), | Dimerization, a 
“al. /¢ al fs s < 
cal./gm. cal. gm. a __ cal /gm. S$ 
639.3 —3.6 0 é 
a 
604.9 —3.0 —2.6 o 
q 
564.8 —8.1 —7.0 Zz 
- - c = ° 
$225 —15.3 —13.7 4 
fr 
178.1 —24.6 —22.0 
4 
$33.1 —34.8 —31.8 oO 
= 
385.3 —47.8 —42.3 
333.9 —64.2 —54.2 
283.7 —78.9 —68.1 
231.6 —94,9 —86.0 a 
Fi 
171.8 —114.3 —110.7 
se 
ol 
struction of the enthalpy-concentration chart are pre- Ww 
sented as supplementary material in Table 3 and Figure in 
6. At 760 mm.Hg. the system forms a maximum boiling ef 
azeotrope at 77.4 wt.c¢ formic acid and 107.2°C. in ill 
agreement with the data of Minjer and Reinders reported 4 
fi 


by Horsley ‘'®), Two sets of data were obtained as formic 
acid tends to decompose when mixtures of high acid Fe 
content are subjected to prolonged refluxing. The first 


PABLE 3 


Vapor-Liguip EQUILIBRIA DATA FOR THE SYSTEM ForMIC 
\cip-WATER AT 760 MM. HG. 


7 ° ; 
Temperature, | x, Mole Fraction 
c. Acid in 
Liquid Phase — | 


y, Mole Fraction 
Acid in 
Vapor Phase 





Refluxing Time: 2 hr. 
101.0 0.047 0.020 
101.6 0.098 0.047 
102.1 0.156 0.082 
103.4 0.223 0.124 
104.1 0.281 0.174 
105.5 0.370 0.268 
106.0 0.412 0.332 
107.1 0.567 0.552 
107.2 0.573 0.573 
106.4 0.674 0.733 
105.0 0.792 0.853 
104.8" 0.794 0.847 
103.8 0.843 0.908 
102.7 0.901 0.952 
101.7 0.948 0.978 
Refluxing Time: 1 hr. 
105.2 0.768 0.839 
104.8* | 0.794 0.851 
103.5% 0.857 0.920 
103.0 0.882 0.944 
102.5 0.913 0.962 
100.9 0.962 0.989 





*Estimated by graphical interpolation 





1.0 TABLE 4 
THE EFFECT OF REFLUXING TIME ON ForMIc AcID 














09 DECOMPOSITION 
— @ REFLUXING TIME = 
Q 08 © REFLUXING TIME : aims Refluxing Liquid Vapor 
nae = | jp Time, | Composition, Composition, 
o O7 hours wt. % acid vt. % acid 
o = —_ —_ = - - - - ———————————— 
i, a 0.5 90.8 93.6 
3 0.6 ; 
a ee 1.0 90.8 93.6 
oq 08 2.0 90.7 93.4 
4 
Bs 04 3.0 90.8 93.3 
o - - 
5.0 90.6 93.0 
Q 0.3 
& 7.0 90.0 92.6 
On Bae 
an 
$ 
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LETTERS TO THE EDITOR 





Menlo Park, California. 
Editor, C.J.Ch.E.: 


RE: The Effect of the Less Viscous Liquid in the 
Laminar Flow of Two Immiscible Liquids by 
T. W. F. Russell and M. E. Charles published in 
the February 1959 issue of C.J.Ch.E. 


This paper brings to mind some unpublished work 
which L. R. Handley and the writer did in 1952 and 
which culminated in U.S. Patent 2,821,205 of January 
1958. 

At the time, we were considering means of reducing 
the pressure drop of pipe lines pumping extremely highly 
viscous crude oil. By means of injecting a thin film of 
water at the pipe wall, we were able to accomplish the 
effect which is so well illustrated in the paper by Messrs. 
Russell and Charles. 


In order to minimize the amounts of water to be 
handled and the chance that this water might emulsify 
with oil, we restricted ourselves to amounts of water less 
than 2% of the crude oil pumped. It is interesting to note 
that in a laboratory installation with over 50 feet of 
approximately 1-inch pipe there appeared to be no mixing 
of oil and water and substantially the entire water film 
could be removed intact at the terminal of the line. 


Two problems became apparent for the practical appli- 
cation of this method in long crude oil pipe lines. The 
first is the necessity for cleaning the pipe so thoroughly 
that all traces of oil are removed from the pipe wall and 
then treating this wall in such a way that it will remain 
water-wet. If this is done successfully, even temporary 
stoppages of the water flow will not cause oil to wet the 
wall, and the water film can be re-established with ease. 
But except for new clean lines, the process necessary to 
remove all traces of oil is extremely tedious and expensive. 


The second problem is concerned with the difference 
in specific grav ity between oil and water. Fortunately, 
this difference is seldom great, and highly viscous crude 
oils tend to be closer to the gravity of water than less 
viscous ones. Nevertheless, in long lines there will be a 
tendency to thin the annular layer of water on one side 
of the pipe, to thicken it on the other, and thereby 
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reduce the efficiency of the process. If the difference in 
the specific gravity of the fluids is very small, that of the 
water could be economically adjusted by the admixture 
of other fluids. 


It may be of interest to note that if the water film is 
very thin compared to the diameter of the pipe and 
when the viscosity of the oil is very large compared to 
that of water one can visualize the oil as flowi ing like a 
solid plug lubricated by the water film in the pipe. More- 
over, it is a fair assumption in this case that the velocity 
in the water film varies nearly linearly from zero at the 
wall to that of the plug at the interface. In this case the 
thickness of the water film divided by the pipe diameter 
will be approximately 3 the ratio of the quantity of water 
pumped to that of oil. This can be verified by div iding 
Equation 22 by Equation 21 in the paper and using the 
assumptions made above. 


E. G. Chilton, 
Stanford Research Institute. 


Edmonton, Alta. 
Editor, C.J.Ch.E.: 


We thank Dr. Chilton for his comments and agree that 
7 Qw r,—r 
the relation — 
Qo m 
the original paper) is approximately valid for plug flow 
under the conditions he describes. 


(using the nomenclature of 


Dr. Chilton mentions that the wetting characteristics 
of the pipeline wall and the density differential between 
oil and water phases are two of the problems to be con- 
sidered in the commercial application of the concentric 
flow system. Current pilot pipeline investigations at the 
Research Council of Alberta with oil and water phases 
of equal density indicate that for a given oil there is a 
region of stable concentric flow provided the pipe is 
maintained in a water-wet condition. 


M. E. Charles and T. W. F. Russell, 
Research Council of Albei ‘a. 
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INSTRUCTIONS TO AUTHORS 


Manuscript Requirements 


i. 
z 


Tables should be numbered in 


The manuscript should be in English or French. 


The original and two copies of the manuscript should 
be supplied. These are to be on 83 x 11 inch sheets, 
typewritten, and double spaced. Each page should be 
numbered. 


Symbols should conform to American Standards 
Association. An abridged set of acceptable symbols is 
found in the third edition of Perry’s Chemical E ngi- 
neers’ Handbook. Greek letters and subscripts and 
superscripts should be carefully made. 


Abstracts of not more than 200 words in English indi- 
cating the scope of the work and the principal findings 
should accompany all technical papers. 


References should be listed in the order in which they 
occur in the paper, a text, using the form 
shown here: “Othmer, F., Jacobs, Jr, J. J., and 
Levy, J. F., Ind. Eng. vt 34, 286 (1942). Abbrevi- 
ations of journal names should conform to the “List 
of Periodicals Abstracted by Chemical Abstracts”. 
Abbreviations of the common journals are to be 
found, in Perry’s Handbook also. All references should 
be carefully checked with the original article. 


Arabic numerals. They 
should have brief descriptive titles and should be 
appended to the paper. Column headings should be 
brief. Tables should contain a minimum of descriptive 
material. 


All figures should be numbered from 1 up, in Arabic 
numerals. Drawings should be carefully made with 
India ink on white drawing paper or tracing linen. All 
lines should be of sufficient thickness to reproduce 
well, especially if the figure is to be reduced. Letters 
and numerals should be carefully and neatly made, 
preferably with a stencil. Generally speaking i originals 
should not be more than twice the size of the desired 


reproduction; final engravings being 3} in. or 7 in. 
wide depending on whether one column or two is 
used. 


8. Photographs should be made on glossy paper with 
strong contrasts. Photographs or groups of photo- 
graphs should not be larger than three times the size of 
the desired reproduction. 


9. All tables and figures should be referred to in the text. 


Submission of Manuscripts 
1. The three copies of the manuscript, including figures 
and tables, should be sent directly to: 
DR. A. CHOLETTE, editor, 
The Canadian Journal of Chemical Engineering, 
Faculty of Science, Laval University, 
Boulevard de |’Entente, 
Quebec, Que. 


The authors addresses and titles should be submitted 
with the manuscript. 


The author may suggest names of reviewers for his 
article, but the selection of the reviewers will be the 
responsibility of the editor. Each paper or article is to 
be reviewed by two chemical engineers familiar with 
the topic. Reviewers may request that they remain 
anonymous. 


All correspondence reviews should be 


directed to the editor. 


regarding 


Reprints 
1. At least 50 free 
supplied. 


“tear sheets” of each paper will be 
Additional reprints may be purchased at cost. An 
estimated cost of reprints, with an attached order 
form, will be sent to the author with the galley proofs. 
Orders for reprints must be made before the paper has 
appeared in the Journal. 
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